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CHAPTER  I 


INTRODUCTION 


Fracture  mechanics,  the  study  of  failure  mechanisms 
and  the  development  of  failure  prediction  criteria  for  en- 
gineering structures,  focuses  considerable  attention  on 
the  stress  distribution  in  the  neighborhood  of  a crack  em- 
bedded in  an  elastic  body.  The  intensity  of  the  elastic 
stress  field  around  the  tip  of  a sharp  crack  plays  a domi- 
nant role  in  fracture  analysis.  Irwin  (1)  showed  that  the 
Griffith  (2)  fracture  theory,  originally  based  on  energy 
considerations,  is  equivalent  to  an  approach  using  the  con- 
cept of  crack-tip,  stress  intensity  factors.  The  basic 
idea  of  the  combined  theory  is  that  the  crack  propagates 
under  static  load  when  the  stress  intensity  in  the  region 
near  the  crack  tip  reaches  a certain  critical  value,  say 
Kc«  For  a given  material  the  combination  of  crack  size 
and  loading  at  incipient  fracture  results  in  a stress  in- 
tensity factor  equal  to  Kc«  When  the  load  or  crack  size 
is  kept  below  this  threshold  point,  the  severity  of  the 
crack  tip  is  measured  by  the  "stress  intensity  factor"  K 
which  assumes  the  general  form 

K = 0f  *Y*  VTT  & 

This  K factor  is  a quantity  determined  analytically  or 
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experimentally  and  varies  as  a function  of  the  crack  length 
a,  the  crack  configuration  and  external  geometry  factor  y , 
and  the  manner  in  which  the  internal  loads  or  far  field 
stresses  a ^ are  applied.  Successful  application  of  the 
theory  requires  that  the  crack  tip  plastic  zone  size  be 
small  compared  to  the  crack  dimensions.  This  permits  the 
use  of  elasticity  theory  to  determine  the  stress  intensity 
factors . 

As  a rule,  the  crack-tip  stress  state  in  a structural 
component  must  be  described  by  at  least  two  parameters , say 
Kj  and  Kjj.  Kj  corresponds  to  Mode  I loading  which  pro- 
duces symmetric  stress  states  near  the  crack  tip  and 
stands  for  the  amplitude  of  the  skew- symmetric  portion  or 
Mode  II  loading.  The  criterion  of  fracture  now  requires  a 
combination  of  Kj  and  reaching  some  critical  value,  i.e. 


F (KJf  Ki;r) 


F 

c 


If  the  stress  state  is  triaxial  in  nature,  a third  parameter 
KIIJ;  is  usually  required  and  the  fracture  criterion  must 
be  modified  accordingly. 

Common  experience  indicates  that  a crack  follows  a 
path  in  a structural  component  as  dictated  by  the  stress 
field.  At  each  instant,  the  direction  of  crack  propagation 
depends  on  the  energy  state  and  the  material  properties  in 
a region  ahead  of  the  crack  tip.  It  is  of  utmost  importance 


2 


to  be  able  to  predict  the  direction  of  crack  initiation 


One  common  method  used  to  predict  crack  growth  which 
may  be  useful  in  mixed  mode  problems  is  the  criterion  of 
maximum  strain  energy  release  rate.  The  strain  energy 
release  rate^  represents  the  amount  of  energy  released 
during  an  incremental  increase  in  crack  area.  The  critical 


strain  energy  release  rate  criterion  and  the  Kc  criterion 


can-  be  related.  For  specimens  subject  to  Mode  I loading, 

» j k ^ 

;£j  equals  C.  I where  C.  is  a constant  depending  on  whether 
E 1 

the  body  is  in  a state  of  plane  stress  or  plane  strain. 


Analytical  problems  for  the  stress  intensity  factors 
are  available  for  a number  of  crack  geometries  and  boundary 
configurations  (3-11) . The  largest  class  of  solved  prob- 
lems apply  to  the  conditions  of  Mode  1 and  II  loading  for 
states  of  plane  stress  or  plane  strain.  A large  complement 
of  mathematical  methods  have  been  used  to  seek  solutions 
and  are  reviewed  in  (9,10).  The  major  analysis  techniques 
used  for  K calculations  are  complex  variable,  conformal 
mapping.  Green's  function,  Riemann-Hilbert  problem,  Wiener- 
Hopf  method,  singular  integral  equations,  and  asymptotic 
approximations.  Authors  utilizing  these  methods  have  pro- 
duced hundred?  of  solutions  many  of  which  are  tabulated  in 
(11) . Analytical  solutions  for  three-dimensional  cracked 
bodies  are  more  limited,  generally  restricted  to  the  penny- 
shaped and  elliptical  crack  geometry. 

Analytical  solutions  for  K factors,  although  large 
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in  number,  do  not  fully  satisfy  the  designer's  needs  when 
analyzing  complex  structures.  The  presence  of  stiffeners, 
irregular  geometries  and  complex  boundary  loadings  often 
complicate  the  real  problem.  Thus  considerable  attention 
has  been  focused  on  numerical  methods  particularly  finite 
element  techniques. 

Early  studies  of  fracture  mechanics  problems  involv- 
ing the  finite  element  method  have  been  carried  out  by 
Swedlow(12),  Tuba  (13),  and  Watwood  (14).  They  attempted 
a straightforward  application  of  the  technique  with  no 
special  attention  given  to  the  stress  singularity.  The 
specimen  was  conceptually  divided  into  finite  elements  with 
a relatively  high  eler^nt  concentration  near  the  crack  tip. 
Computational  experiments  indicated  the  need  to  include  a 
very  large  number  of  degrees  of  freedom  in  order  to  obtain 
reasonable  accuracy  (12-14) . The  numerical  value  of  the 
stress  intensity  factor,  determined  in  this  manner,  varies 
over  a considerable  range,  depending  upon  which  node  is 
chosen  for  its  calculation.  Further  study  of  this  technique 
(15,16)  concerns  questions  about  its  convergence,  reliabil- 
ity, and  convenience. 

In  an  attempt  to  eliminate  some  of  these  undesirable 
features,  an  alternative  approach  has  been  developed  by 
Wilson  and  others  (15)  and  extended  by  Oglesby  and  Lomacky 
(16) . They  directly  incorporated  both  the  finite  element 
method  and  the  analytical  crack  tip  expansions  to  form  a 


special  singular  crack  tip  element.  This  element  centered 
at  the  crack  tip,  has  the  deformation  field  given  by  the 
elastic  singular  solution,  with  the  stress  intensity  factors 
and  crack  tip  displacement  components  as  the  parameters 
characterizing  its  behavior.  The  crack  tip  element  is  con- 
nected to  standard  elements  along  its  boundary,  by  requir- 
ing the  nodal  displacement  components  on  this  interface  to 
be  consistent  with  the  crack  tip  solution.  Finite  element 
results  for  complicated  geometries  using  special  cracked 
elements  are  believed  to  easily  give  accuracy  within  2-3% 
(17)  . 

Evaluation  of  stress-intensity  factors  will,  in  most 
cases  require  a numerical  attack  due  to  the  complicating 
factor  of  geometry.  The  boundary  integral  equation  (B.I.E.) 
technique  is  a numerical  method  particularly  suited  to  the 
stress  analysis  of  simply  or  multiply  connected  homogeneous, 
elastic  bodies. 

The  method  as  developed  by  Rizzo (18)  exploits  the 
strong  analogy  between  potential  theory  and  classical 
elasticity  theory.  Unlike  Jaswon(19)  Rizzo's  approach  is 
not  limited  to  two-dimensions  and  works  directly  with  as- 
signed boundary  quantities  instead  of  stress  functions. 

This  approach  depends  on  the  knowledge  of  the  singular  sol- 
ution to  the  Navier-Cauchy  elasticity  equations  in  two- 
dimensions  corresponding  to  a concentrated  force.  The 
singular  solution  gives  rise  to  a vector  identity  similar 
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to  Green's  third  identity  for  Laplace's  equation.  Taking 
the  field  point  to  lie  on  the  boundary  contour,  a boundary 
formula  is  obtained  which  gives  a relation  between  boundary 
displacements  and  corresponding  boundary  tractions.  Since 
either  of  these  boundary  quantities,  in  principle,  deter- 
mines the  other,  the  formula  provides  a constraint  between 
them  which  generates  a set  of  simultaneous  integral  equa- 
tions involving  the  unknown  boundary  data.  The  reduction 
to  a problem  involving  only  the  boundary  surfaces  reduces 
the  dimensionality  by  one  and  has  no  requirement  for  dis- 
cretization of  the  region  contained  within  the  boundary 
surface.  This  reduction  in  dimensionality  significantly 
reduces  the  computational  requirements. 

This  method  has  been  extensively  developed  by  Rizzo, 
Shippy,  Cruse  and  Lee  (20-28).  They  have  applied  the  B.I.E. 
method  to  the  following  class  of  problems;  1)  elasto- 
dynainics , 2)  three-dimensional  elasto-static,  3)  elastic 
inclusion,  and  4)  anisotropic  elastic  boundary  value  prob- 
lems . At  the  present  time  application  to  cracked  bodies 
is  restricted  to  symmetric  geometries  which  allow  the  prob- 
lem to  be  changed  from  a cracked  body  to  a mixed  boundary 
value  problem  with  no  cracks.  The  unsymmetrical  crack 
problem  is  excluded  in  Rizzo's  formulation  due  to  regular- 
ity requirements  on  the  boundary  surfaces.  These  regular- 
ity conditions  require  that  the  body  be  finite,  and  smooth, 
and  representable  by  unique  parametric  expressions. 


Classical  fracture  problems  require  that  the  crack  faces  be 
defined  to  lie  on  the  same  contour  and  meet  sharply  at  the 
crack  tip.  These  conditions  violate  the  regularity  require- 
ments upon  which  Rizzo's  formulation  relies. 

It  is  the  purpose  of  this  thesis  to  extend  the  Rizzo 
technique  to  permit  its  application  to  those  crack  problems 
lacking  symmetry.  Chapter  2 presents  the  mathematical  for- 
mulation of  the  B.I.E.  for  plane  cracked  bodies,  while  the 
solution  technique  is  applied  to  a specific  Mode  I,  Mode  II 
problem  in  Chapter  3.  Chapter  4 presents  results  and  dis- 
cussion for  two  sets  of  problems.  The  first  set  is  for 
verification  of  the  solution  method  and  the  second  presents 
data  for  the  bent  crack.  Finally,  Chapter  5 is  concerned 
with  conclusions  and  applications  of  the  developed  method. 


CHAPTER  2 


MATHEMATICAL  FORMULATION 


Solutions  of  boundary  value  problems  in  fracture 
mechanics  require  the  analysis  of  bodies  containing  mathe- 
matical models  of  crack-like  flaws.  The  geometric  modeling 
of  these  flaws  generally  result  in  configurations  or  geo- 
metries that  are  non-regular  as  defined  in  (29,30).  The 
basic  integral  equation  formulation  developed  by  Rizzo  (18) 
requires  that  the  body  under  consideration  be  finite  and 
bounded  by  a single  smooth  contour  C , the  contour  satis- 
fying certain  regularity  conditions  (29)  . One  regularity 
condition  requires  that  contour  C possess  a unique  rep- 
resentation in  the  parametric  form  x = x (s).  A sec- 

a a 

ond  condition  is  imposed  on  the  surface  derivatives  x^(s) 
which  must  be  continuous.  The  two  crack  faces  in  classical 
fracture  problems  are  defined  to  lie  on  the  same  contour; 
thus,  unique  parametric  representation  of  the  individual 
surfaces  is  excluded.  The  presence  of  a reentrant  vertex 
at  the  crack  tip,  and  the  existence  of  sharp  corners, 
where  the  crack  intersects  the  outer  boundary,  violates 
the  second  regularity  condition.  Section  2.1  presents  a 
summary  of  the  Rizzo  formulation  considering  C to  be 
regular.  This  permits  solution  of  a class  of  problems  ex- 
cluding the  elastostatic  crack  problems  arising  in  fracture 


mechanics.  The  mathematical  processes  required  to  overcome 
these  restrictions  and  to  extend  the  B.I.E.  method  for 
solution  of  general  crack  problems  are  the  subject  of 
Section  2.2 

2 . 1 Basic  Integral  Equation  Formulation 

The  integral  equation  formulation  presented  in  this 

section  is  restricted  to  the  analysis  of  classical  plane 

elastostatic  problems  for  which  the  material  may  be  taken 

1 2 

to  be  isotropic  and  homogeneous  ' . The  Navier  displace- 
ment equations  of  equilibrium  in  the  absence  of  body  forces 
are  given  by 


a + ^VafS  +)i  Vcm  ' 0 (“'S  * 1'2> 

(2.1) 

where  the  displacement  is  denoted  by  u^  (x^)  while  Xfl 
are  the  orthogonal  cartesian  coordinates  x^>x2  an<*  ** 
are  the  familiar  Lame  constants.  Eq.  (2.1)  formulated  for 
problems  of  plane  strain,  may  also  embrace  generalized 
plane  stress  with  suitable  exchange  of  material  constants. 
The  solution  to  this  differential  equation  must  also  satis- 
fy the  boundary  conditions  for  the  displacements  and 


■^A  more  general  formulation  for  the  3-D  elastic  prob- 
lem may  be  found  in  References  21,25,27. 

2 

The  formulation  presented  in  this  section  is  very 
similar  to  that  of  Rizzo(18).  It  is  presented  here  because 
it  is  the  foundation  for  mathematical  extensions  germane 
to  this  thesis. 


r 


tractions  on  the  contour  C , given  by 


u$  lxa>  " qB 


\ <=  O' 


and 


a Qn  = tQ  (x  ) * e C 

a$  a 8 a a 


(2.2) 


1 2 

where  C , C are  those  portions  of  the  contour  subject  to 


displacement-traction  loading  boundary  conditions  respect- 
ively. The  unit  vector  n is  the  outward  normal  vector 


for  the  body  D with  components  nQ . 


The  body  under  consideration  will  be  assumed  finite 
and  bounded  by  a single  smooth  contour  C,  which,  as 
described  in  (29) , admits  a unique  representation  in  the 


parametric  form  x = x (s) . Further,  the  derivatives 

a a 


x'a(s)  are  assumed  to  be  continuous.  The  parameter  S 
is  the  arc  length  along  the  contour  from  some  arbitrary 
origin. 

Figure  1 illustrates  a simply  connected  body  with 
closed  contour  C.  For  the  body  (D  + C)  under  the  action 


of  boundary  tractions  t , the  relation 


= X UB,8  nc+  unB  <uc,B+  “B.c1  (2.3) 


must  hold  on  C. 

The  distance  between  field  points  Q with  coordinates 
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rtltir'lfilhrf  md*  II  t v f I II  - mi 


R 


I 


xa  and  the  source  points  p with  coordinates  iQ  is  given  by 


r = [ (x  - £ ) (x  -£  ) ] 

a a a a 


(2.4) 


Then,  concentrated  forces  of  magnitude,  -47ru(X  +2y)/(X  +3y)  , 

are  applied  at  points  p = p(?a)  in  the  xa  directions  res- 

* 

pectively,  giving  rise  to  displacement  vector  u~  as 

p 

defined  by  Love  (31) 


where 


U3  “ Ua6  Sa  = U18  + U2B 


UaS  f ^aB^  r + M r'a  r'3 


(2.5) 


and  where  ea  are  unit  base  vectors.  Eq.  (2.5)  satisfies 
Eqs.  (2.1)  in  D,  provided 


r ^ 0,  M = 


(X  + u) 

X+  3y 


(2.6) 


In  Eq.  (2.5)  and  in  what  follows,  all  differentiation 
is  with  respect  to  the  field  point  x.  That  is. 


- (x  -?  ) 
r a 


3x  na 


r (x0  'Va 


where  the  normal  is  evaluated  at  xa  also.  The  traction 

* * 
vectors  tg  on  C corresponding  to  up  are  computed 

from  Eq.(2.3)  and  may  be  written 


The  vectorial  form  of  Betti's  reciprocal  work 
theorem  may  be  written 

/(u  T Q - t U Q)dS  = 0 (2.8) 

a aB  a aB 

C+m 

where  dS  is  an  element  of  arc  length  at  Q and  where  the 
point  p = P(£a)  of  D has  been  excluded  by  a small  circle 
m of  radius  p because  of  the  singular  nature  of  UQ  and 
T 0 for  r = 0.  The  vectors  and  their  derivatives  are 
taken  to  be  non-singular  and  sufficiently  continuous  for 
the  usual  validity  of  the  theorem.  These  vectors  corres- 
pond to  an  equilibrium  stress  state;  body  forces  associated 
with  U q and  T 0 

~p 


M » 


are  taken  as  zero.  It  can  be  shown 


(18)  that 


Lim 
p-*-o  m 


Ac  Uci 


dS  = 0 , Lim 

p-*o 


fu  T 
J a a 

m 


0 as  = -«  "1»8(p) 


(2.9) 


where 


a = (X+3u)/(4ttu(X+2w))  . 


Therefore,  Eq.  (2.8)  yields 

ue  (p)  = a / [ua  (Q)  T^ip.O)  - ta  (Q)  uoe(p,0)  ]dS(Q) 
c 

(2.10) 

which  is  the  plane  counterpart  of  Somigliana's  identity  (31) 
for  the  displacements  inside  the  body,  D,  caused  by  surface 
tractions  and  displacements.  Eq.  (2.10)  is  analogous  to 
Green's  third  identity  of  potential  theory  which  expresses 
a harmonic  function  in  terms  of  the  boundary  values  of  the 
function  and  its  normal  derivative.  A parallel  develop- 
ment using  Green's  third  identity  leads  to  the  formulation 
of  the  Mode  III  problem. 

The  properties  of  the  integrals  in  Eq.  (2.10)  as  the 
boundary  C is  approached  are  necessary  for  the  development 
of  B.I.E.  relating  the  displacements  and  tractions  on  C. 

The  integrals  present  in  Eq.  (2.10)  are 

<Mp)  = ft  (Q)  Ua8(p,Q)  dS'  VP*  = /ua  (Q)  TctS(P'Q)dS 

p c M c 
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Lim  <J>g(p)  = ipg(P)  = it*  fK-2p  M]Ug  + f ua  ^ Tag*P,Q^dS 
P P ° (2.11) 

This  result  is  readily  verified  provided  ua  is  assumed  to 
satisfy  a Holder  condition  on  C and  provided  the  singular 
integral  is  evaluated  in  the  sense  of  a Cauchy  Principal 
Value  (32) . A similar  limiting  procedure  for  the  first  in- 
tegral gives 

Lim  <J>ft(p)H  <J>fl  (P)  = ft  (Q)  U r(P,Q)  dS 
p+P  3 3 5 “ aB 

(2.12) 

The  integral  equations  relating  the  surface  tractions 
to  surface  displacements  are  found  from  Eq.(2.10),  to  (2.12). 
Since 


U6(P)  = (P)  ~ VP)] 


(2.13) 


and 


air  (K  - 2yM)  = 1/2 


it  results  that 


uB  (P)  - 2 a /ua  (Q)  T B(P,Q)dS(Q) 
p c p 

’ *2  2 /*« 


(Q)  Uag(P,Q)dS(Q) 


(2.14) 
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As  shown  by  Rizzo  (18)}  Eq.  (2.14)  will  be  regular 
according  to  the  classification  of  singular  operators  by 
Muskhelishvili  (32)  for  values  of  Poisson's  ratio  us [0,1/2). 
Given  tfl  on  C,  the  traction  problem,  the  two  components  of  Eq. 
(2.14)  form  a simultaneous  pair  of  integral  equations  for 
the  unknown  u^.  The  rigid  body  displacement  vector 


UB  <*>  " aS  + beae3  5a<P> 


satisfies  the  homogeneous  form  of  these  equations,  i.e. , 
ug(P)-2  S J“ Ja  (Q)  Tag(P,Q)dS(Q)=  0 (2.15^ 


c 

where  a^  and  b are  arbitrary  constants.  Application  of  the 
Fredholm  alternative  yields 


L 

J a 


/G 


(Q)  dS(Q)  = yE3Baxg(Q)  ta(Q)  dS(Q)  = 0 
c 

(2.16) 


as  the  necessary  and  sufficient  conditions  for  the  solu- 
bility of  Eq.  (2.14).  These  are,  of  course,  the  conditions 
that  the  body  be  in  equilibrium. 

It  should  be  noted  that  although  the  Fredholm  alter- 
native was  employed  above,  Eq.  (2.14)  is  not  a Fredholm 
integral  equation,  since  the  term 


T*g  = K (Ln  r),a  nr(Ln  r),^ 


is  present  in  the  kernel  function  T 


a6 


More  precisely. 


15 


i 


#. 





whereas  a limit  of  zero  is  required  for  Fredholm  kernels. 

Nevertheless,  the  index  k in  Muskhelishvili  Chapter  19  (32) 

computed  from  T 0 equals  zero  and  the  Fredholm  alternative 
exp 

is  applicable. 

Given  the  u^  on  C,  the  displacement  problem,  the  two 
components  of  Eq.  (2.14)  form  a simultaneous  pair  of  sin- 
gular integral  equations  for  the  unknown  t on  C.  Rizzo 
proved  in  (18)  that  the  homogeneous  equations 

/ta(Q)  Uag(P,Q)dS(Q)  = 0 
C 

has  no  nonzero  solutions.  Thus  it  may  be  expected  that  the 

displacement  problem  formulation  of  Eq.  (2.14)  results  in 

unique  solution  ta  on  C for  arbitrary  IT  on  C. 

With  the  knowledge  of  u.  * uJ^  on  part  of  C and 

P P 

(2) 

fca  = t„  on  the  remaining  part  (mixed  problem) , a set  of 
four  equations  on  the  unknowns  u^^  and  t^^  on  C is  ob- 
tained. The  mixed-mixed  boundary  value  problem  may  also 
be  formulated  by  extension  of  the  above  reasoning.  Evalu- 
ation of  the  displacement  field  (and  subsequently  the  stress 
and  strain  fields)  internal  to  the  boundary  C may  be  ob- 
tained through  Eq.  (2.10)  from  the  boundary  data  given  on 


all  contours. 


2 . 2 Integral  Equation  Formulation  for  Cracked  Plate 
Geometries 

In  this  section,  the  Boundary  Integral  Equation  (B.I.E.) 

formulation  necessary  to  solve  elastostatic,  edge  cracked 

3 

plane  problems  is  developed  . This  formulation  is  restric- 
ted to  traction  loading , although  displacement  or  mixed 
boundary  conditions  may  be  addressed  with  equal  ease. 

Direct  formulation  of  the  Boundary  Integral  Equation 
(z.14)  applied  to  edge  crack  geometry  leads  to  an  incomplete 
problem  statement.  This  shortcoming  is  caused  by  the  inher- 
ent restriction  that  the  complete  contour  of  body  D be  des- 
cribed by  a unique  parametric  representation  x = xQ  (s).  The 

p p 

upper  and  lower  crack  faces  violate  this  condition.  The 

other  restriction  on  the  continuity  of  x'  (s)  must  also  be 

p 

addressed  since  a discontinuous  change  in  (s)  occurs  at 
the  crack  tip,  and  may  occur  at  the  crack  corners. 

The  method  is  first  illustrated  for  the  case  of  a 
crack  with  Mode  III  loading,  which  is  equivalent  to  the 
Dirichlet  problem.  This  classical  problem  (30,33)  is  formu- 
lated beginning  with  Green's  third  identity  yielding  the 
governing  integral  equation 

W(P)  - £ /w(0)  * K(Q,p)dS(Q)  = F (p) 

C (2.17) 

3"I  Throughout  this  thesis  attention  is  restricted  to  edge 
crack  problems.  However,  the  method  discussed  is  appli- 
cable to  internal  crack  geometries  by  reformulating  Eq.(2.14) 
to  account  for  multiply  connected  regions. 


i 


34K 


where : 


"F(P)  ~hj Ik  <“>“  ro  as<°> 

C 

k(0,p)=  (x-yfw*  (y-y0>^ 

r2 

o 

■r 

r o - i <x'xo) 2 + (y-y0> 2i 

Eq.  (2.17)  relates  the  unknown  displacement  W in  the  z- 
direction  to  an  applied  shear  load  which  is  proportional  to 
^ (Q)  on  the  boundary  C.  Q and  P are  field  and  reference 
points  on  the  boundary  C,  and  n is  the  unit  normal  at  point 
Q,  and  x,y  are  the  usual  cartesian  coordinates.  The  deri- 
vation of  Eq.  (2.17)  requires  the  same  restrictions  as  im- 
posed on  Eq.  (2.14)  pertaining  to  unique  parametric  repre- 
sentation of  C and  continuous  surface  derivatives  x^(s)  • 
Because  the  Mode  III  governing  integral  Eq.  (2.17)  is 
simple  and  has  similar  mathematical  behavior  to  the  Mode  I 
and  II  formulation  it  is  selected  for  initial  investigation. 

A simple  Mode  III  B.V.P.  is  defined  with  geometry 
illustrated  in  Fig.  2.  It  consists  of  a right-circular 
cylinder  of  radius  R with  a sharp  2^^  wedge  cut  out  begin- 
ning at  the  center  line  of  the  bar.  The  selection  of  an 
appropriate  traction  loading  results  in  a closed  form  ana- 
lytical solution  suitable  for  comparison  with  the  integral 
equation  results.  Analytical  solution  of  the  Mode  III 
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problem  begins  with  the  definition  of  the  usual  cylindrical 
coordinates  r,0,z  with  origin  at  the  bar's  center.  The  gov- 
erning differential  equation  is  the  classical 

V2  W = 0 

where  W = W ( p , 0 ) p = r/_ 

/ R 

0<P<1  , il)  <0< 

The  boundary  conditions  are  defined  as. 


a)  on  the  upper  wedge  face,  0 =ipQ 


y 1 3y  = T =0 
P 9 0 z 

b)  on  the  lower  wedge  face,  9 = 2 tt  -i/r 


c) 


on 


u I »L  = T = 

v P 30  z 


the  circular  boundary, 


3w 
3 P 


= T = ^ cos 
P =1  z R 


P = 1 

( ) 
( 2^"V  j 


Solution  of  the  harmonic  equation  and  satisfaction  of  the 
B.C.,  yields  the  following  results  for  the  boundary  displace- 
ments 


a) 

on 

0 = 

lb 

yo 

b) 

on 

0 = 

2tt-4> 

c) 

on 

P : 

= 1 

W(p  ) = PB 
t o 

w(p,-4jq)  = -pB 
W(1,0)  = cos  J B ( 0 -ipQ) 


where 


B = tt/2  (tt  -ipQ) 
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Integral  Eq.  (2.17)  is  formulated  for  a modified  cir- 
cular wedge  geometry.  The  crack  tip  and  crack  corners  are 
rounded  with  circular  arcs  of  radius  to  insure  starting 
with  a contour  C whose  surface  derivatives  are  continuous 

and  ip  is  assumed  to  be  nonzero  so  the  two  crack  faces  are 
o 

separate  and  distinct.  This  modified  geometry  shown  in 
Fig.  3 meets  all  the  restrictions  imposed  on  Eq. (2.17) . 

The  upper  wedge  surface  is  labeled  C^,  the  circular 
boundary  C2,  the  lower  wedge  surface  and  the  crack  tip 
and  upper  and  lower  corner  contours  , Cj-  amd  Cg  respec- 
tively. Eq.  (2.17)  is  applied  to  this  modified  geometry 
and  the  resulting  equation  separated  into  six  parts.  Each 
of  the  parts  represents  a unique  portion  of  the  I.E.  for 
reference  points  located  on  individual  contours  , 

i = 1 to  6.  This  yields: 


Wi(Pi)  " ? / W1(Q1)K1(Q1,Pi)dS.(Q1) 

C1 

“if  /W2(Q2>K2<Q2Pi>dS2(Q2>  - ? /v03>K3(03'Pi>dS3<°3) 
c 2 C 3 

-?  /« 4<°4)K4(Q4,Pi,dS4(04)-  f J” 5<°5)K5<°5,Pi>dS5<05> 
C4  c5 

- 1 A 

TT  J € 


■6(WQ6'Pi)dSfi<06> 


Fi(Pi)  i = 1,6 


no  sum  on  1 


(2.18 
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where : 


:if  dS3<Q3> 


Kj  (Qj’Pi>  = 


<x3  * xoi>|si  + <yj  - yoilfj4 


ij  " {<xj  ' xoi>2  + (yj  ' yoi)2  } 


w^  = displacements  on 


P^,Q^  = field  and  reference  points  on 


Points  and  include  all  points  in  the  set  which  form 
the  contour  C^;  dS^  are  incremental  arcs  on  Ci  for  i=l  to  6. 
The  operation  limit -* 0 is  performed  on  Eq.  (2.18).  The 
integral  over  the  surfaces  C 4,  and  Cg  when  i = 1,2,  or  3 
yields: 


Lim  / W4(Q4)K4(Q4,P.)dS4(Q4)  = Lim  ^W-  (Q;.  ,p.)  dS5  (Qs> 
C4  2 C5 

/wK6<°  's'Vas/V  0 i = 1-2.3 


dS6<06)  ‘ 0 i - 2-2-3 


provided  points  P at  the  crack  tip  and  each  of  the  cor- 

K 

ners  are  excluded  from  the  set  of  reference  points  P^.  The 


k 


exclusion  of  points  from  Eq.(2.18)  imposes  no  limitation 
on  the  eventual  numerical  solution  of  the  equation  set.  This 
is  an  obvious  result  of  the  numerical  solution  technique 
which  requires  discretization  of  the  equations  with  the  ac- 
companying freedom  of  selecting  a finite  number  of  reference 
points.  The  elimination  of  a small  set  of  reference  points 
p offers  no  special  restriction  nor  any  great  inconvenience. 
The  exclusion  of  points  P also  allows  the  remaining  integral 
equations  for  i = 4,5  and  6 in  the  Eq.(2.18)  system  to  be 
ignored.  For,  when  the  limit  er-0  is  performed  and  the  points 

p are  excluded  from  the  reference  set,  these  equations  are 
K 

eliminated . 


I 


-A 


sin  (9n  ~ ip0) 


(A)  to'  d<?  + w /ft  \ 

$2  -2^* cos  (0 -A)+l  W2  ^®o 


2ir'*o  1 

- I / w2(e,  -id8  + 1/w 


sin (9q  + yQ)  d6 


3 6 ^~26  *cos  (6  +i|»  )+l 

o o 


= T7i^0)  /cos°-^C^fS 

0 

1 / ,,,  50  si] 
if  J wi(S)  rr-ns  s0 


. Ln  I 2-2cos (9 


-v 


(2.19b) 


6o  sin  2 <l>0 


1 /* 

TF  / W.(9)  jrl^-cos(Q^Q)  2 d9  + W,(6  ) 
y 2'  l-26n  cos (9+to0) +6-/  3 o 


2tt-i|> 


- 4Ti%/c°=  2^fcf  L" 


•cos  (9+i^q)  +6{ 


(2.19c) 

In  Eq.(2.19)  V^,  V*2,  and  W3  are  the  unknown  displacements 
in  the  z-direction  for  the  surfaces  C^,C2,  and  C3  respec- 
tively, and  6,6q,9,9o  are  the  remaining  cylindrical  co- 
ordinate variables  for  the  field  and  reference  points.  The 

points  6 =1,  6 = 0,9  = ii/  , and  9 = 2n-  \1>  are  inadmis- 

o o o o o o 

sible  reference  points  in  the  Eq.  (2.19),  and  ip#  e(0,7r)  . 

Numerical  solutions  of  the  integral  equations  which 
arise  in  this  paper  require  particular  attention  to  those 
singular  kernels  whose  integrals  exist  only  in  the  sense  of 
a Principal  Value.  Special  techniques  to  account  for  the 


1 m 
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effect  of  the  singular  stress  field  near  crack  tips  and 
other  critical  points  are  also  necessary  for  good  numerical 
accuracy.  Therefore,  a major  section  in  Chapter  3 is 
devoted  to  the  discussion  of  these  techniques,  and  numeri- 
cal results  pertinent  to  this  chapter  are  presented  without 
explanation  of  the  numerical  processes  used  to  obtain  them. 

Eq.  (2.19)  is  numerically  solved  for  selected  values 
of  0 £ £ 90°.  Displacement  values  alonq  the  wedge  faces 

are  of  particular  interest.  A convenient  parameter  helpful 
to  the  evaluation  of  the  numerical  results  is  obtained  by 
dividing  the  computed  displacements  along  the  wedge  faces 
by  the  exact  solution,  and  fitting  a least  squares  straight 
line  to  this  data.  The  intercept  value  for  5-  0 is  labeled 
NDP  (Normalized  Displacement  Parameter) . With  the  use  of 
the  parameter  NDP,  all  problem  solutions  for  0 <_  ipQ  <90° 
will  yield  NDP  equal  to  one  if  the  numerical  error  is  zero. 
Percent  deviation  of  NDP  from  unity  is  thus  a means  of 
evaluating  the  numerical  accuracy  of  the  results . The  re- 
sults are  tabulated  in  Table  2.1. 


TABLE  2.1 

Anti-Plane  Wedge  Problem  Results 


*o 

m 

BS 

20° 



15° 



| 

B 

NDP 

btj 

.999 

1.02 

BE 

D 

i 

i 

SB 

.2% 



.1% 

.i% 

.6% 

2% 

10% 

-v20% 

HBSfll 

f 


The  results  obtained  for  ipQ>  30°  clearly  support  the 
conclusion  that  the  I.E.  (2.19)  is  a proper  representation 
of  the  sharp  wedge  problem  and  may  be  solved  accurately. 

For  values  of  fQ<  30°,  the  important  displacement  results 
near  the  crack  tip  develop  significant  errors  as  ap- 
proaches zero.  These  results  are  so  badly  behaved  that  they 
render  NDP  meaningless  when  »!>-*■  0 and,  the  determination  of 
the  stress  intensity  factors  for  = 0 is  impossible  by 
this  formulation. 

The  requirement  for  a unique  parametric  representation 
of  all  boundary  surfaces  causes  the  computational  difficul- 
ties experienced  above.  A direct  limiting  process  executed 
on  Eq.  (2.19a)  and  Eq.  (2.19c),  taking  special  care  with  all 
Principal  Value  integrals,  shows 

, Lim  (Eq. 2.19a)  = . Lim  (Eq.  2.19c) 

^o  -*-0  vo  -*■  0 

This  result  is  clearly  manifested  in  the  numerical  example. 
The  determinate  of  the  linear  algebraic  equation  set  ap- 
proaches zero  as  approaches  zero,  confirming  that  Eq. 
(2.19)  is  incomplete  in  the  limit  il>  + 0 . 

Before  proceeding  with  the  developments  necessary  to 
complete  the  I.E.  it  is  convenient  to  rewrite  Eq.  (2.19)  in 
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source  points  on  the  upper  and  lower  wedge  faces  are  cal- 
culated from  Eqs.  (2.20a)  and  (2.20c).  As  0-*O  the  co- 
efficients become  equal,  one  approaching  the  limit  value 
from  above  and  the  other  from  below.  This  suggests  that 
if  a sum  and  difference  state  relating  adjoining  wedge 
face  source  points  is  formed  then  higher  order  functions 
necessary  to  complete  the  governing  I.E.  may  be  identified. 
A scalar  6q  which  measures  the  distance  from  the  crack  tip 
to  the  paired  source  points  is  introduced.  The  sum  and 
difference  state  is  achieved  by  first  adding  and  then  sub- 
tracting Eqs. (2.20a)  and  (2.20c)  from  each  other  with  the 
result : 


(«!<«„>  + V6o>)  *i / W|KlllslV  + K13<VSo)|dSl 


-1/ 


w2(s2) 


K21 *S2'6o)+K23  *s2'6o^ 


ds. 


W3  *s3* 


K11<s3'5o)+  Ki3<s3'V 


ds. 


F1<V  + F3(V 


W2*S2}  " iAcSlf^l^)  dsl  “it  f' 

ri  r, 

Si' 


(2.21a) 


W2 (s2) K22 (s2' s2^dS2 


W3(s3)K32(s3'S2)  dS2  = F2(s2} 


(2.21b) 
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UIMB* 


(W1<V  - W3<So>)  - i/W  1<S1> 


KU<s1'So>-K13 


<*l'4o>]  ds. 


-l/w 


2 (S2} 


K21<S2'V  + K23<S2'V  ds2 


r2 


- f/V'af1 


K31 ^S3,5o^“  K33 


r3 


(s3f6)  j 


ds. 


- Fl(5o}  - F3<V 


(2.21c) 


The  operation  Lim  t (j  -»  0 performed  on  Eqs . (2.21a)  and 
(2.19a)  leads  to  identical  equations.  Taking  the  limit 
as  4*q -*  0 in  Eq.  (2.21c)  results  in  the  null  identity.  This 
limiting  character  of  Eq.  (2.21c)  suggests  the  following 
assumption  concerning  the  dependence  of  the  integrals  upon 
^ 0 . Assume, 


?/Wl<«l>[Kll<»l'So>'  K13<sl'V]dSl  - Wl<{<,H*onfl«o>] 

n>l  as  'J'q+O 

;/W2<s2>[K2  1<«2'4o>'  K23(s2'4o)]dS2*°|  *>"£2<V  I 


n > 1 as  <|M-0 


?/W3<=3)[K31<=3'40»-  K33's3'4o>jdS3'-M3(61>)+0|  +of3f4»)) 


Fl(4o'*o)-  WV  ’ 0 


n 

4>n  F 


(8.)]  "i 


n > 1 as  4)^0 


1 as  ^0  (2.22) 


3 

.jb 

I 
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Then  Eq.  (2.21c)  assumes  the  form, 


W - W - / W+  C fi(5o>j 

+ C • W + (M3<V  **no  ■ £3<5o>Wno  • f(V 

' (2.23) 

and  after  simplifying  (2.23)  it  follows  that 

C + £2<V  + W - 5(So>]  * 0 

(2.24) 

Eq.  (2.24)  may  be  satisfied  if  equals  zero  or  by  the 

bracketed  quantity  vanishing. 

In  the  following  it  is  assumed  that  the  bracketed 
expression  equals  zero,  thus  implying  the  existence  of  an 
additional  constraint  equation. 

The  arguments  presented  in  Eqs.(2.22)  to  (2.24)  do 
not  provide  a convenient  method  for  extracting  the  neces- 
sary equation,  although  its  existence  is  implied.  Assuming 
n is  integer,  it  is  more  convenient  to  differentiate  n 
times  and  take  the  limit  as  thus  arriving  at  the 

bracketed  quantity  directly.  Application  of  the  technique 
shows  that  only  one  differentiation  is  required  for  the 
problem  considered  here. 

Returning  to  the  notation  of  Eq.  (2.19)  this  leads 
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to  the  simple  process 


L 


Lim  ^ l / Lim  rv 

— ( Eq.  (2.19a)-Eq.  (2.19c)}  = ^0^  [Eq.2.21c] 


V03+o  ) 


(2.26) 


Executing  these  operations  yields 


Lim 


3W,  ,6 


V°  ( Ho 


^'U0/w  i(S» 


6 *sin  2^o 
o 


6 2— 2 6 * 6 *cos2i|;  +6  2 
o o o 


d6 


3W3(6o) 

3^ 


cos(0-ii>  ) 1-6  ‘cos  ( Q ) 

o o o 

cosTo-lFT^ffi- 2 ~ 1-26  'Cos(0-i|>  )+6  2| 
o Yo  o o ro  o 1 


d0 


*h-  (' 


w3(S) 


60  sin  2<p 


dfi 


2t T-lp( 


6 -26*6  *cos  2 * +6 

O o o 


.±.2_  f 

2 J 


Tz  Ln< 


<J>o 


i-26  *cos (0-tf>  ) +6 


1-26 q* cos (0+^Q) + 6q* 


)d0 


(2.27) 


Eq.  (2.27)  will  be  shown  to  provide  the  additional  infor- 
mation necessary  to  complete  the  I.E.  representation  of 
the  closed  crack  anti  plane  problem.  In  the  following 
a numerical-analytical  technique  is  developed  to  allow 
refomulation  of  Eq.  (2.27)  into  a non-singular  integral 
equation  thereby  permitting  direct  numerical  solution. 

The  technique  begins  by  identifying  those  integrals  which 
exist  only  in  the  principal  value  sense  in  the  limit  as 

0 ; in  this  case, the  integrals  over  the  crack  faces 
containing  and  The  Cauchy  singular  kernel  in 
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the  integral  precludes  interchanging  the  limit,  differ- 
entiation, and  integral  operators.  To  proceed,  the  sing- 
ular points  in  each  P.V.  integral  are  identified  and  the 
integrals  separated  into  three  integration  domains  as 
follows : 


P.V.  f W (5) 
°6  + H 

+P.V.  J Wx  (6) 

6-h 


6-h 


K (6,t|»0)  d6 


-A 


(6) 


1 

K(6,^q)  d6  +/wi(6> 

6+h 


K(6,V<0)  d6 
K(5^0)  d6 


(2.28) 


For  example,  K (6  may  have  the  form 

. <5  *sin  2i|» 

K(5/^0)  a o 12 

62-26 • 6 -cos  2^  +6  2 
o To  o 

K (6  ,\p  ) is  singular  for  6=  6 in  the  Lim  \h  0.  The 
o o ro 

integrals 

6-H  1 

/Wi(6)  K(6,Viq)  d6  , y*W1(6)  K (6,1^)  d« 

0 6 + h 

are  Riemann  integrable  and  sufficiently  well-behaved  to 
permit  interchange  of  differentiation  and  limit  operations 
with  the  integral.  The  P.V.  integral  in  Eq.(2.28)  has 
been  restricted  to  the  domain  of  integration 


r 


6q  - H <_  6 <_  6q  + H . The  value  H is  selected  to  be  real 

and  positive  such  that  it  is  always  less  than  or  equal  to 

the  smallest  value  of  6 . Since  the  numerical  solution  of 

o 

the  integral  equations  proceeds  by  the  method  of  finite 

differences  requiring  satisfaction  of  the  equations  at  a 

finite  number  of  collocated  points,  H is  conveniently 

chosen  as  the  value  of  the  division  size. 

* 

/ Further  analysis  of  the  P.V.  integral  requires  an 
approximation  of  the  functional  behavior  of  the  displace- 
ments on  the  crack  surfaces.  The  & functional  depend- 
ence near  the  crack  tip  is  well  documented  (34,35).  A 
two  term  approximation  for  the  crack  tip  displacements  in 

the  domain  6 - H <6  < 6 + H is  selected  which  generally 

o — — o 

takes  the  form: 


= C1  <5  + C2 


For  surfaces  near  the  crack  tip,  a value  of Y = 1/2 
is  used.  The  linear  formy  = 1 is  used  for  surfaces  removed 
from  the  crack  tip. 

The  constants  and  C2  are  dependent  on  the  values 
of  y selected  as  well  as  the  values  of  the  unknown  dis- 
placements evaluated  at  points  6 = - H,  and  + H. 

The  Cauchy  principal  value  integral  may  now  be 
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approximated  by 

- 60  +H 


V.  (6,<J>Q)  d6  “ Cx  P.V.Jly' 


6 o ~H 


fio+H 


+ c2*p 


• V.  y”  K(6,i^o)  d6  - 


6o-H 


K (6,4»o)d6 


6n-H 


cr  W + C2*J2<V 


The  P.V.  integrals  J.^  j and  J2(ip  ) are  analytically 

evaluated.  Differentiation  w.r.t.  followed  by  Lim  0 

o 

are  operations  easily  performed  on  the  resulting  closed 

form  expressions  for  ci  Ji  ft  ) + C2  J2 ft  ) * Analytical 

expressions  result  which  depend  on  the  value  of  Y and 

are  called  insert  functions.  Performing  the  remaining 

operations  in  Eq.(2.27)  yields  the  desired  equation. 

A system  of  integral  equations  which  result  from 

Lim  (Eq.2.21)  may  now  be  written.  Using  the  detailed 
'J'o-O 

notation  of  Eq.  (2.19)  and  completing  the  details  for  Eqs. 
(2.21a,b)  yields 


- wi<V  + ? f w; 


1-6  *cos  0 

o 

1-26  *cos  6+6  ‘ 
o o 


f Tz  (G>*  Ln  j l-26o*cos0+^l  2d0, 


-W3(6o) 


(2.29a) 
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=)!>-•  =*!*-■ 


»in  9, 


COS0  +1 
o 


d0  " w2(0o) 


z*^11  r cos 

+ k J V”  - kj  w3  I 

o o *■ 

= TZ  f2  Tz(0)  Ln  j2-2  cos  <0“0o)j  d0 


cosG  +1 
o 


(2.29b) 


from  Eq.  (2.21c)  it  follows  that 

0 _+H 
O 


, 2/5“  / \ |<Z5“-/5“^ 

W (6)  i - !■—  (/6+H  + /6  -h)+  Ln — 

1 77  H \ ° ° / (/5“+/5“^ 

o o 

. /*27r  (1-6  2)  Sin  9 

I / w (9)  —5 2 — d0 

71  J * (1-26  COS0  + 5 2) 

o o o ' 


-/S^H)  (/5^+/6o+h> 


/S^H)  (ZS^H-ZS? 


26  COS0+6  2) 
o o 


+ W3(6o} 


. 2 /5“  / _ \ ( 

- — (/67+h  + /6  -H)+  Ln  - 
TT  H \ O ° / j 


(/T^-ZT^H)  (ZT>  /6q+H) 

(ZT>ZT^h)  (Z6q+h-ZS^) 


o -h  _ i r 1 

/°  26  . f 26 

W,(6)  ^ <36  " J I W (6)  j 

u*-‘-,aJ  y+„  bs-4o>  J 

O 

26^  Sin  0 

2 d9 

1-26  CosG  + 6 
o o 

m 


(2.29c) 
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In  the  foregoing 

me  [0,1]  9 e 10 , 2tt] 

Pe  [0 , 1]  B o e [0 , 2ir] 

and  Tz  are  normalized  tractions  prescribed  on  the  cir- 
cular boundary.  Eq.  (2.29c)  incorporates  the  insert 
function  assuming  piecewise  square  root  displacements  on 
the  crack  surface,  i.e.y  = j . 

Investigation  of  the  set  of  integral  equations  (2.29) 
shows  them  to  be  a quasi-regular  system  of  Fredholm  equa- 
tions of  the  second  kind  as  defined  in  Muskhelishivili , 
pp.  415-417  (32).  This  assures  that  a unique  solution 
exists.  It  remains  to  show  by  example  that  this  solution 
represents  the  cracked  anti-plane  boundary  value  problem 
of  elasticity. 

A numerical  solution  of  Eq.  (2.29)  was  obtained.  The 
details  of  this  numerical  technique  shall  be  discussed  in 
a special  section  in  Chapter  3.  Comparison  with  the  right 
circular  cylinder  (4»  * 0)  problem  discussed  earlier,  for 

loading  tz(m)  = \ Cos  j tractions,  provides  the  neces- 
sary verification.  The  numerically  determined  displacement 
field  on  the  boundary  surfaces  were  accurate  within  1%  as 
compared  to  the  analytical  solution.  These  results  con- 
firm the  limiting  process  developed  for  the  anti-plane 
crack  problem.  Eq.  (2.29)  may  now  be  assumed  to  be  a 
complete  formulation  of  the  anti-plane  cracked  body 
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problem  of  elasticity. 

2.3  The  Plane  Elastostatic  Clacked  Body  Integral  Equations 
The  plane  elastostatic  cracked  body  integral  equations 
are  developed  by  extension  and  generalization  of  the  tech- 
niques used  in  the  anti-plane  problem.  A traction  loaded 
plane  body  with  m number  of  edge  cracks  as  shown  in  Fig. 4a 
is  considered.  Fig.  4b  illustrates  the  modified  "open" 
crack  geometry  necessary  to  begin  the  analysis.  This  geo- 
metry retains  the  sharp  crack  tips  and  corner  point 
description.  Each  open  crack  surface  and  all  external 
boundaries  are  modeled  by  continuous  parametric  expressions 
which  best  lend  themselves  to  the  shape  under  consideration. 

A single  parameter  ‘I'q  is  chosen  which  allows  modeling  the 

Xiiin 

modified  geometry  of  Fig.  4b  and,  for  the  ^ — g'  y^elds  tlie 

closed  crack  geometry  of  Fig.  4a.  The  Kth  crack  surfaces 

are  labeled  C^,  where  the  upper  (u)  surface  has  been 

arbitrarily  chosen  to  differentiate  it  from  the  opposite 

crack  face,  called  the  lower  (L)  face.  The  set  of  all 

crack  surfaces  C?  + cj'  + • • • + + • •+  cj*  + cj* 

11  K K mm 

is  designated  C.  The  external  boundary  is  labeled  S. 

The  plane  elastostatic  system  of  I.E.  (2.14)  applies 
directly  to  the  modified  geometry  with  the  exclusion  of 
crack  tip  and  corner  reference  points.  The  exclusion  of 
these  points  eliminates  the  need  for  rounding  the  corners 
which  introduces  additional  complexities.  As  in  the  anti- 
plane problem  no  theoretical  compromise  or  numerical 
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inconvenience  occurs  from  this  restriction. 

This  coupled  pair  of  I.E.  is  separated  into  parts  by 
writing  individual  equations  for  each  crack  surface  and 
external  boundary.  Recalling  Eqs . (2.5)  and  (2.7)  for 


definitions  of  and  Eq.  (2.14)  assumes  the  form 


- 2 -f  TaMPK'  Qr)  dS  ^ 


u 


- 2 


i/"c.  (Ok)  as<0K) 

'L. 


- 2 


“ KcQU)  Tae  <p“. 


Q)dS(Qu) 


^ u _U  « U . 

C1  + C2  + CK-1+  CK+1+"*CM 


- 2 


“ /“a<0L)  TcS^k'O)  dS(o1') 


C^+C^t  +cu  +cu  + cu 

12  K-l  K+l  * M 


-2  “ y*ua(Qs)  Tag(PK'Q>  ds<Qs> 

- JVQ)  Ua6(PK,Q)  dS(Q) 


"S 

= - 2 


C+S 
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(2.30P  uK) 


rum  n i -m  &'**'*«& 


W - 2 o/v<»  Tce(ps'Q>  ds<°> 

C 

-2  ayu0(as)  TclB(Ps,as)  ds(es) 

5 

- - 2 e.A(a)  uae(Ps.a)  as(o> 


<V  T«e<Ps'as>  dS<«s> 


= - 2 


ue (pk»  - 2 


Jo.e<ps'a>  ds(a> 


-h 


<°K>  Tc,B<PK'QK)  dS(QK> 


2.  K(QK>  To6  <PK' 


q£)  dS(Q^) 


(2.30eS) 


a /u  (QU)  T 

-J  ° 


'.6(PK'  Q)  ds(a“) 


cu+cu+  +CU  +CU  + CU 

12 * ’ * K-l  LK+1  ’ ’ *Cm 


a /v aL)  J6  (^K-a)  <js  (qli 

Cl+C2+' ' * +CK-1+CK+1+  *+Cm 


“-/“a 


(as)Ta8(PK'0)  dS(Qs> 


-2  i /\,<a>UaB(P J,  0)  ds  (Q) 


(2.30  gLK) 
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where 


a = 1,2  3 = 1,2  K = l,m  (no  sum  on  K) 

Eqs.  (2,30)  may  be  considered  as  a system  of  coupled 
integral  equations.  The  K index  represents  the  decomposi- 
tion of  the  I.E.  over  the  selected  crack  faces.  Eqs. 
(2.303UK)  represent  the  subset  of  equations  whose  reference 
points  lie  on  upper  crack  surfaces,  there  being  K = l,2..m 
such  equations,  one  for  each  crack.  In  like  manner,  Eq. 

(2. 303LK)  represents  m equations  whose  reference  points 
lie  on  the  lower  surfaces.  Eq.  (2.303S)  is  the  remaining 
equation  for  reference  points  located  on  the  external 
boundary  S.  As  an  example  of  this  notation,  the  function 
Up(p£)  is  interpreted  as  the  displacements  in  the  x^ 
direction  for  the  set  of  points  which  lie  on  the  upper 
surface  of  the  Kth  crack.  The  kernel,  T -(P^,Q^)  is  the 
function  defined  in  Eq.  (2.7)  evaluated  for  reference  points 
on  the  upper  Kth  crack  surface  (P^)  and  field  points  on  the 

i\ 

lower  Kth  crack  surface  (Q^) . Integration  variables  con- 
taining sub  or  superscripts  assume  their  normal  meaning. 

The  rewriting  of  Eq.(2.14)  as  Eq.  (2.30)  permits 
the  separation  of  the  integrals  into  Cauchy  principal 
value  and  Riemann  integrals.  In  both  Eq.  (2.303uK)  and 
Eq.  (2.303LK),  the  P.V.  integrals  appear  as  the  first  two 
integrals  in  the  equation,  the  remaining  integrals  being 
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Riemann  integrals.  In  Eq.  (2.306S),  the  integrals  over  the 
boundary  S exist  only  in  the  P.V.  sense,  all  others  are 
Riemann  integrable. 

The  similarity  of  Eq.  (2.30)  to  the  governing  I.E. 
for  the  anti-plane  problem  is  evident.  It  may  be  shown 
that  the  kernel  functions  in  Eq.  (2.30)  and  those  in  the 
anti-plane  representation,  Eq.  (2.19),  have  similar  be- 
havior. The  kernels  are  still  Cauchy  singular  although  the 
complexity  of  the  kernel  functions  increases  greatly. 

This  result  could  be  anticipated  since  both  formulations 
rest  upon  the  same  mathematical  development.  Following 
a similar  approach  as  in  the  Mode  III  problem  a complete 
set  of  governing  I.E.  is  formulated  for  Mode  I and  II  prob- 
lems. Proceeding  as  in  the  Mode  III  case,  Eq.  (2.30)  is 
recombined  and  a operation  performed  in  the  following 

sequence. 


^ iXg  { Eq.  (2.240UK)  + Eq . ( 2 . 246LK) J 
° (2.31a) 

^ J1®  { Eq.  (2.24BS)}  (2.31b) 


Lim  - , 

'J'q-  ° { Eq.  (2.240UK) 


Eq. ( 2 . 24$LK)> 

' (2.31c) 


where  6=1,2,  K=  1,2.  . .m 

The  above  operations  are  performed  for  each  6,  K 
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equation  with  special  care  exercised  for  those  integrals 
which  have  Cauchy  singularities.  Eq.  (2.30a)  and  (2.30c) 
each  yield  a minimum  of  2m  equations  and  Eq.  (2.30b)  yields 
a minimum  of  two  more  equations.  Thus,  a geometry  with  m 
external  cracks  requires  a minimum  system  of  4m  + 2 coupled 
integral  equations.  Complicated  crack  or  external  boundary 
shapes,  requiring  multiple  parametric  representation  may 
require  further  decomposition  of  the  equations  with  a com- 
mensurate increase  in  the  system  of  governing  equations. 

An  example  problem  is  solved  using  the  geometry  of 
Fig.  2 for  \pQ  = 0 and  for  surface  tractions  resulting  in 
both  Mode  I and  Mode  II  states.  The  detailed  results  are 
reported  in  the  next  chapter  and  compare  within  1%  to  the 
exact  analytical  solution.  The  accuracy  of  the  solution 
of  Eq.  (2.29)  for  the  anti-plane  Mode  III  example  and 
Eq.  (2.31)  for  the  plane,  Mode  I,  II  case  gives  a numerical 
verification  of  the  effectiveness  of  the  method. 
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CHAPTER  3 


MODE  I AND  II  PROBLEM  STATEMENT  AND  INTEGRAL  REPRESEN1 


The  solution  technique  for  Mode  I,  II  and  III  load- 
ing developed  in  Chapter  2 requires  numerical  verification. 
A circular  disk  with  a simple  edge  crack  is  selected  for 
study  because  it  possesses  a known  closed  form  analytical 
solution  for  certain  specific  boundary  conditions.  This 
disk  is  in  a state  of  plane  strain  and  subjected  to  Mode  I 
& II  loading. 

An  accurate  solution  of  this  problem  via  Eqs.  (2.31) 
serves  as  a numerical  verification  of  the  plane  elastostat- 
ic  solution  method.  A bent  crack  geometry  which  has  drawn 
considerable  attention  [36,37]  recently  for  application  to 
mixed  mode  fracture  is  also  considered  to  demonstrate  the 
general  applicability  of  the  method.  This  geometry  is  asy- 
metric  and  involves  stress  singularities  at  the  bend  vertex 
in  addition  to  those  at  the  crack  tip.  A limiting  case  of 
this  geometry  is  the  straiqht  crack;  therefore  only  the 
general  bent  crack  formulation  needs  to  be  presented,  since 
both  straight  crack  and  bent  crack  results  follow.  The 
bent  crack  problem  statement  and  equation  formulation  is 
presented  in  Sec.  3.1,  and  the  development  of  the  numerical 
solution  technique  is  discussed  in  Sec.  3.2.  The  mode  III 
problem  was  extensively  discussed  in  Sec.  2.2.  It's  numer- 
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ical  solution  is  obtained  by  a similar  process  to  that  pre- 
sented in  Sec.  3.2  and  need  not  be  discussed  further. 

3.1  Bent  Crack  Problem  Formulation 

A cracked  circular  domain  of  radius  R,  as  shown  in 
Fig.  5,  is  considered.  A cartesian  coordinate  system,  xy, 
having  its  origin  at  the  disk  center, is  employed.  A 
straight  crack  of  length  d intersects  the  circular  surface 
and  lies  on  the  x axis  in  the  region  xe[f,R].  At  x=f  a 
secondary  straight  crack  of  length  c emanates  at  an  angle  a 
to  the  x-axis.  The  values  of  c,  d,  and  a are  restricted 
only  by  the  requirement  that  the  secondary  crack  must  not 
penetrate  the  circular  boundary  S.  Solutions  for  the  trac- 
tion boundary  value  problem  for  several  loadings  on  the  ex- 
ternal surfaces  are  to  be  obtained  numerically. 

Analysis  of  this  problem  begins  with  the  construction 
of  the  modified  open  crack  geometry  shown  in  fig.  6.  The 
secondary  crack  is  opened  from  the  crack  tip  by  an  angle 
2 v|j0  so  the  primary  crack  has  uniform  separation 
2g=2£* cosa* tan  The  x-coordinate  value  of  the  crack  tip 

is  labeled  e while  the  projection  of  c on  x is  denoted  by 
£„  The  individual  surfaces  forming  the  boundary  are  iden- 
tified by  T.  i=l,2***5  (see  Fig.  6).  A normalized  geom- 

1 r 

etry.  Fig.  7,  is  obtained  by  dividing  all  length  variables 
by  the  radius  R,  capital  letters  denoting  the  normalized 
variables.  On  the  bent  crack  surfaces  r,  r_  the  distance 

X / D t 

from  the  crack  tip,  p,  locates  field  points,  while  B iden- 


[ 

L 


tifies  source  point  locations.  On  the  straight  crack  sur- 
faces T.  field  points  are  located  by  the  variable  xl 

and  source  points  by  B2 , both  measured  from  the  coordinate 
origin  . On  the  circular  boundary  r_  the  radian  mea- 

J / 

sure  0 locates  field  points  and  0Q  the  source  points.  As 
previously  discussed  the  modified  geometry  needs  no  crack 
tip  or  comer  rounding  if  the  set  of  reference  variables 


are  restricted  to  an  otien  set  excludincr  these  points. 

The  pair  of  coupled  B.I.E.  (2.14)  apply  directly  to 

the  modified  bent  crack  geometry  with  crack  tip  and  corner 

points  excluded  from  the  reference  set.  Eqs.  (2.14)  are 

separated  into  five  parts  by  writing  individual  equations 

for  reference  points  located  on  the  discrete  I\  boundaries. 

This  decomposition  permits  examination  of  the  P.V.  integrals 

and  allows  application  of  the  solution  method  discussed  in 

Chapter  2.  Thus,  Eq.  (2.14)  yields. 

C/  . 

— cos  ip 

-1  r 0 

-j-UjfP)  + /[u1(p)  • KIj^P,p)  + vx(p)  • KIIj;L]dp 

0 

+ y[u2(xl)  KIj2(P/3tL)  + v2(xl)  • KIIJ2]  d xl 

F+C*  tan  ip  • tana 
o 

27t-C*tan 

. o 


I>*|H 


J 4 * KIJ4<P,51)  + v4(xl)  • KIIJ4  ]dx} 


F-C*tan  ip  • tan  a 
o 


C/ 

— cos  ip 


+/tu5(p)  • H„<P.»  + v5,P)  . J 


dp 


Fj  (P) 


C/ 

— cos  \p 


(3. 1AJ) 


VJ  (P)+  j{VP>  ■ KXHJl(p,p)  + Vi< 


P)  • KIV 


J1 


|dp 


X 

yju2  (STL).  KIIIj2(p,x1)  + v2  (xl)  • KIVJ2  | d icl 


F + C • tan  ^Q*tana 


2-rr-Otan  tpQ 

/{u3  (S)-  KIIIJ3(P,9o)  +v  (0)  . KIVj3  l d0 

C tan  90  ^ 


Hu,  <*1>.  KIIIj4  <P,*1)  + v<  (xl)  . KIVj4  n 
F-C*tan  ip  • tan  a 


u 


c/ 

— ' cos  ip 


FtB  (P) 


KIjl(P,Q)e{k  - 4.m.M.RJL  x 

L 


(3.1BJ) 


KIIjl  (P,Q)  = |-4«yM*RJL 


• \ 8LnRJL 

y JL'xf  ~9 JT~ 


K ‘ { (Ln  rjl> 'y  V(Ln  rjl) ,x  ny ^ 

, . t \ 9LnR 

Kill  (P,Q)  5 J-4*  )j*M*  Rtt  • Rtt  1 — 2L- 

JL  ( JL,y  JL,x  f 9nT 


K|[(Ln  Rjl)  ,y  nx  - (Ln  Rjl)  ,x  yy]| 


Kivjl  (p'Q)  = |[K-4*  p* 


W-PTT  Rtt  ] \ 
JL,y  JL ,y J f 


9 Ln  R. 


RJL  - j(XJ  (Q)-  XL(P))2  + (YJ(Q)  - VP))2  p 
nT  = unit  normal  on  surface  r_ 

Li  Li 

M = - (X+y)/(X+3y) 

K = 2y2/(-X+3y) 

X 5 -(X+3y)/2iry(X+2y)  ) 

y,X  are  Lame's  material  constants. 

In  the  notation  of  Eq.  (3.1  AJ)  and  (3.1BJ),  J takes 
on  values  from  1 to  5 creating  a system  of  ten  coupled 
equations.  The  reference  point  variable  P is  understood  to 
have  meaning  B,  B£  or  0q  depending  on  the  reference  sur- 
face for  which  the  equation  applies.  The  kernel  functions 
KIjl,  Kill jk i KHIjl'  an<*  kIVjl  defined  in  Eq.  (3.1),  are 
presented  in  algebraic  detail  in  Appendix  A.  Kernels  hav- 
ing repeated  indices  in  Eq.  (3.1  AJ)  and  (3.1  BJ)  are  singu- 
lar in  the  lim  \p0  -+  0 and  integrals  containing  such  kernels 
exist  only  in  the  Principal  Value  sense. 

To  proceed,  form  a sum  and  difference  state  on  those 
equations  having  source  points  on  the  crack  surfaces.  This 
is  followed  by  the  limiting  operation  lim  -+Q,  after  dif- 
ferentiation with  respect  to  in  the  case  of  the  differen- 


lj 

i 

f 

t 

~ 
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ce  equations.  These  operations  may  be  summarized  as  fol- 
lows : 


\ 


L 


Y " 0 (Eq.  (3.1  A 1)  + Eq.  (3.1  A 5)} 

0 

, Lim  (Eq.  (3.1  A 2)  + Eq.  (3.1  A 4)} 

^ Lim0  {Eq.  (3.1  A 3) } 

Llm  {-jra-lEq.  <3.1  A 2)  - Eq.  (3.1  A 4)]} 
V»0  - O dV0 

^oLimo  {a|^tEq*  (3-1  A 1)  - Eq.  (3.1  A 5)]} 

* LJ“0  {Eq.  (3.1  B 1)  + Eq.  (3.1  B 5)} 

Lima  {Eq.  (3.1  B 2)  + Eq.  (3.1  B 4)} 

. L^m  {Eq.  (3.1  B 3) } 
r0 

ib  Limo  (3* 1 B 2)  - Eq.  (3.1  B 4)  ] } 

^oLimo  3^{lEq-  <3-1,B  - Ec3 • (3.1  B 5)]} 


(3.2  A 1) 

(3.2  A 2) 
(3.2  A 3) 
(3.2  A 4) 

(3.2  A 5) 
(3.2  B 1) 
(3.2  B 2) 
(3.2  B 3) 
(3.2  B 4) 
(3.2  B 5) 


The  mathematical  operations  indicated  in  equation 
(3.2)  are  performed  using  the  techniques  discussed  in  Sec- 
tion 2.2.  The  singularities  are  removed  by  assuming  a 
piece-wise  functional  behavior  for  the  unknown  displace- 
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ments  followed  by  analytical  evaluation  in  the  neighborhood 
of  the  singular  point  of  each  principal  value  integral.  A 
key  requirement  for  the  successful  performance  of  this  ana- 
lytical step  is  knowledge  of  the  correct  functional  behav- 
ior of  the  unknown  displacements.  Fortunately,  this  behav- 
ior is  well  known. 

In  general  a satisfactory  representation  is  given  by 
u = C.  ^ +C 2 where  the  selection  of  the  parameter  y de- 
pends upon  the  surface  in  question.  This  selection  of  y is 
treated  in  detail  in  Section  3.2.  The  analytical  evalua- 
tion of  the  integrals  in  the  neighborhood  of  the  singular 
regions  results  in  algebraic  functions  labeled  insert  func- 
tions. The  equations  which  result  are  a coupled  system  of 
simultaneous  non-singular  integral  equations  of  the  second 
kind. 

The  completed  equation  set  (3.2)  is  presented  in  Ap- 
pendix B.  Section  B.l  contains  the  functional  form  of  the 
integral  equations,  detailing  the  geometric  kernel  func- 
tion. Section  B.2  contains  three  sets  of  insert  functions, 
one  for  each  of  three  classes  of  assumed  displacement  func- 
tionals. Section  B.3  presents  the  detailed  forcing  func- 
tions for  general  traction  loading  on  all  surfaces. 

3 . 2 Numerical  Solution  Technique 

The  numerical  solution  of  Eqs.  (3.2)  proceeds  direct- 
ly by  employing  the  standard  methods  developed  for  non- 
singular Fredholm  equations  of  the  second  kind.  A finite 


The  in- 


nuinber  of  mesh  points  partition  the  boundaries  T\ 
tegrals  in  (3.2)  are  similarly  partitioned  and  the  unknown 
displacements  approximated  within  these  intervals.  The  in- 
tegrals are  numerically  evaluated,  resulting  in  a set  of 
influence  coefficients  which  operate  on  the  unknown  dis- 
placement vector  evaluated  at  discrete  mesh  points.  The 
continuous  representation  of  Eqs . (3.2)  is  thus  transformed 
into  a linear  system  of  algebraic  equations. 

The  left  side  of  Eqs.  (3.2)  shown  in  Section  B.l  leads 
to  the  algebraically  equivalent  system  involving  the  influence 
coefficient  matrix  and  displacement  column  vector.  The 
right  side  of  (3.2)  generates  a column  vector  of  forcing 
functions  representing  the  traction  loading. 

[ C ] { u } = { F } (3.3) 

The  resulting  matrix  equation  (3.3)  after  row  elimination  or 
coefficient  substitution  for  the  supression  of  rigid  body 
motion  may  be  solved  directly.  Calculation  of  the  influ- 
ence coefficients  forming  the  [C]  matrix  begins  by  parti- 
tioning the  boundary  as  illustrated  in  Fig.  8.  Ml  points 
are  located  on  and  with  equal  spacing  H1=C/(M1  + .5). 

The  last  reference  point  on  or  r,.  is  located  Hl/2  dis- 
tance from  the  bend  vertex.  Surface^  and  contain  M2  points 
spaced  H2  = D/(M2  + .5)  apart.  The  mesh  point  nearest  the 
bend  vertex  is  again  spaced  H2/2  from  the  bend  point.  A 
total  of  M3  reference  points  are  equally  spaced  on  r 3 with 
H3  = 2tt/^13  + 1)  . No  reference  points  are  located  on  the 


crack  tip,  bend  vertex,  or  crack  corners. 

The  numerical  evaluations  of  the  integrals  proceeds 
by  employing  the  partition  just  described.  The  unknown  dis- 
placements are  approximated  in  a piecewise  manner  by  the 
function 

u = An  CY  + Cn  ye  [Js,l] 

(3.4) 

where  £ is  understood  to  be  6 , xl  or  0 depending  on  the  in- 
tegration domain.  A and  Cn  are  coefficients  involving  the 
unknown  displacements  at  the  mesh  points.  Then,  a typical 
integral  term  in  (3.2)  may  be  written  as 


B-H, 


N 


(n+1)  h. 


/“ i<P)  K <fiP)dp 


n=0 


x(p)  K(p,e)dp 
n*H. 


N 


(n+1)  H 


(n+1) H, 


‘E*-/ 

n=o  n*H. 


prK(p,6)  dp 


+ C"/ 


K(p,6)  dp 

n*H, 


(3.5) 


The  integrals  appearing  in  Eq.  (3.5)  are  further  partition- 
ed into  an  even  number  of  M subintervals  and  numerically 
evaluated  by  Simpson's  rule. 

The  insert  functions  discussed  earlier  and  presented 
in  App.  B.2  contribute  the  dominant  diagonal  coefficients  in 
the  "C"  matrix.  These  functions  result  from  the  analytical 
evaluation  of  the  P.V.  contribution  and  their  accuracy  is 
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directly  dependent  on  the  quality  of  the  assumed  unknown 

displacement  function.  The  non-singular  integrals  in  the 

near  neighborhood  of  the  singular  points  also  result  in 

large  coefficients  located  in  a band  on  either  side  of  the 

diagonal  elements.  The  accuracy  of  these  coefficients  is 

dependent  on  both  the  quality  of  the  assumed  displacement 

function  and  the  accuracy  of  the  numerical  integration. 

The  functional  behavior  of  the  boundary  displacement 

near  the  crack  tip  and  bend  vertex  are  well  documented. 

[36,  38,  39],  The  numerical  analysis  uses  the  square  root 

functional,  y=h,  for  Ml-2  partitions  on  V and  for  Ml-4 

divisions  on  T . The  bend  corner  formed  by  the  common 
s 

point  on  T and  T is  similar  to  a sharp  wedge  139].  The 

4 5 

displacement  field  in  the  neighborhood  of  the  wedge  vertex 
has  the  form  (3.4)  with  y determined  from  the  first  posi- 
tive real  root  of  the  transcendental  equation 

sin  (y[iT+a])  + ysin  (ir+a)  = 0 (3.6) 

Values  of  ae  (0,tr)  result  in  ye  (Sj,l), 

This  displacement  approximation  was  employed  for  numerical 
integration  of  the  non-singular  integrals  and  for  the  ana- 
lytical evaluation  of  the  insert  functions  over  a ten  refer- 
ence point  domain.  This  domain  includes  5 mesh  points  on 
T and  5 points  on  T surrounding  the  vertex  surfaces.  By 

4 5 

restricting  the  bent  crack  angle  a to  be  positive  a linear 
displacement  function  suffices  on  the  upper  bend  vertex. 

This  is  appropriate  because  under  these  conditions  only  the 


lower  surface  possesses  a singular  behavior  near  the  vertex 
[34,39,40].  On  all  of  the  remaining  surfaces,  the  dis- 
placement components  are  modeled  by  piecewise  linear  func- 
tions. It  should  be  noted  that  to  insure  compatibility  of 
results  similar  displacement  functions  were  assumed  for 
both  the  analytically  determined  insert  functions  and  the 
numerically  integrated  coefficients,  depending  on  the  sur- 
face of  integration. 

The  accuracy  of  the  numerical  integration  is  ascer- 
tained by  conducting  numerical  experiments  on  the  most  sen- 
sitive integrals.  Taking  advantage  of  the  16  digit  accur- 
acy of  the  CDC  6600  computer  and  using  the  simple  Simpson's 
formula,  most  integrals  are  successfully  evaluated 
(<<1%  error)  by  a subpartition  of  M = 20.  The  most  sensi- 
tive kernels  require  M = 250.  In  this  manner  the  geometric 
influence  coefficients  forming  [C]  are  calculated.  The 
matrix  [C]  is  full  with  the  dominant  influence  coefficients 
banded  along  the  diagonal. 

Generation  of  the  forcing  function  matrix  proceeds  by 
direct  use  of  the  Simpson  approximation.  Statement  of  the 
traction  loading  on  each  boundary  completes  the  definition 
of  the  integrals  shown  in  App.  B.3  and  evaluation  proceeds 
by  direct  numerical  integration.  A partition  of  two- 
hundred  evenly  spaced  points  for  each  surface  provides 
accurate  results  (<1%  error) . 

54 


.u 


1 


This  numerical  reduction  of  the  integral  equations 
reduces  Eqs.  (3.2)  to  a system  of  2{*2*(M1+M2)  + M3)  alge- 
braic equations  in  the  same  number  of  unknowns.  Solution 
of  this  algebraic  system  by  standard  numerical  methods  then 
leads  to  the  values  of  the  displacement  components  at  the 
mesh  points . 

The  necessary  computations  for  the  examples  presented 
herein  were  performed  on  a CDC  6600  computer.  This  compu- 
ter features  sixteen  digit  accuracy;  thus,  no  special  dou- 
ble precision  operations  were  utilized.  The  program  re- 
quirements dictated  Ml  = 8,  M2  = 7,  and  M3  = 11  as  minimum 
values  for  the  general  bent  crack  problem.  The  upper  limit 
was  established  by  accuracy  requirements  (^1%)  and  seldom 
exceeded  130  equations.  This  equation  size  required 
130,000  memory  units  and  sixty  seconds  to  compute  coeffi- 
cients and  solve  the  equations. 
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CHAPTER  4 

NUMERICAL  EXAMPLES  AND  DISCUSSION 

fl 

Attention  is  now  turned  to  two  classes  of  problems. 

The  first  consists  of  two  problems  for  which  analytical 
solutions  exist,  thus  permitting  verification  of  the 
solution  method.  The  second  class  of  problems  presents 
results  for  the  bent  crack  as  a function  of  crack  angle 
and  loading  mode. 

4 . 1 Numerical  Verification  Examples 

EXAMPLE  #1.  The  circular  disk  with  a straight  edge 
crack  equal  to  the  plate  radius  has  an  analytic  closed 
form  solution  for  specific  traction  loadings.  This  geometry 
is  a special  case  of  the  bent  crack  shown  in  Fig.  5 with 
a equal  to  zero  and  c plus  d equal  to  R.  The  governing 
biharmonic  equation  for  the  Airy  stress  function  is  solved 
using  the  technique  developed  by  Williams  (34,35).  A trac- 
tion loading,  equivalent  to  the  first  term  in  the  Williams' 
solution,  is  thus  generated  to  develop  a special  problem 
possessing  a closed  form  solution.  A problem  summary  of 
the  required  traction  loadings  and  the  resulting  displace- 
ments for  each  boundary  surface  follows  using  the  notation 
in  Chapter  3; 
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u3  (0) 


= _2  kt  l cos6  *sin  (-^-)*(2  (1-v) -sin  9 ) 

2 L l 2 2~ 


-sin  0 -cos  (A)  (cos2  (-|— ) + (l-2v)  ) 


+ 


COS0 


e 

•cos  y~ 


(3 *cos2  |— 


-(2v  + 1)) 


- sin  6 


(3  * sin2|-  + 2v  -4) 


Dimensionless  Displacements  in  the  y-direction: 
vx  (P)  = CQ/p-  KI 

v (xl)=  c Jxl  Ktt 

2 oil 

v4  (XT)  = - co  fsi  kt 

Vc  (p)  = - C Jp  ' K_ 

o o ’ I 

v3  (0)  = ^o  Kxj  sine  sin|-  (2(l-v)-sin2  j-  ) 

+ cos 6 cos  (cos^  + (1-2V)  ) 

^ 1 ) 
(Eq.  continues  on  next  page) 


58 


. 


+ 


sin 


0 


6 2 9 

'cos  — (3  -cos  2~ 


- (2v  + 1)) 


where  C 

o 


4(l-v2) 
E J2~n 


+ 


and 


cosG  • sin  |^.sin2  Jp-  + 2(v-2) 

(4.1) 

v = Poisson's  ratio 


E = Young's  modulus 

and  u,  v represent  displacements  in  the  x,y  directions 
respectively. 

The  loading  for  this  problem  is  confined  to  the  cir- 
cular boundary  and  generates  both  Mode  I and  Mode  II  states 
at  the  crack  tip.  The  stress  intensity  factors  Kj.  and 
are  explicitly  represented  in  both  traction  and  displace- 
ment equations  in  (4.1).  and  K.^  are  arbitrary  and  may 

be  set  equal  to  unity  in  the  traction  equations  above. 

The  B.I.E.  (3.2)  are  solved  for  this  geometry  and 
traction  loading  stated  in  (4.1).  The  resulting  displace- 
ment solutions  are  recorded  and  used  to  calculate  the  stress 
intensity  factors  and  strain  energy  of  the  body.  The  stress 
intensity  factors  are  analytically  defined  by  the  formula 


K 


I 


Lim  Ur,i(  p) 
P + 0 


Co^> 


Lim 

P -*•  0 


Un5( r>  ) 

CQVP 
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'I  " I Will— 


(4.2) 


Lim  Utl.(p  ) _ Lim  Ut5^p  ) 


II  p -*•  0 


P 0 


C /p 
o 


where  Un,Ut  are  normal  and  tangential  displacements. 

The  numerical  calculation  of  these  factors  proceeds 
by  first  forming  the  expressions 


k = Unl  (P)~  Un5(P> 

1(p)  2 C n/P 

o 

77  _ Utl  (p)-ut5  (p) 

2(p>  ‘ 2 cc 

calculated  for  p e rltr5  and  determining  their  values  as 
P + 0 by  extrapolation.  The  calculated  stress  intensity 
factors  then  follow  from  (4.2): 


KI  Kl(0) 


KII  K2  (0 ) 


(4.3) 


The  strain  energy  is  calculated  from  the  boundary  displace- 
ments and  traction  loading  by  use  of  Clapeyron's  theorem(41) , 


- I A »i 


(4.4) 


The  displacements  calculated  at  the  mesh  points  are  numer- 
ically integrated  with  the  tractions  using  Simpson's 
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formula  to  complete  the  calculations. 

EXAMPLE  #1  - Results 

The  results  for  example  #1  are  summarized  in  Table 
4.1  and  4.2.  A parametric  study  is  presented  to  test 
the  sensitivity  of  mesh  size  and  percent  square  root  ap- 
proximation on  the  accuracy  of  the  numerical  solution. 
Fifteen  examples  are  cited  for  fixed  C and  D values  (as  a 
function  of  mesh  size  - Ml,  M2,  and  M3)  and  MT,  the  total 

number  of  equations  in  the  matrix.  The  calculated  values 

_ „ , , , . . Es  calculated 

of  KT,K  T,  and  the  strain  energy  ratio  P— are 

-L  -L-i*  E Exact 

s 

compared  to  unity  to  assess  the  accuracy  of  the  numerical 

solution.  The  maximum  error  in  either  Kj.  , KJ1;  or  Es  is 

listed  in  the  last  column. 

Most  of  the  results  indicate  less  than  1%  error  for 

a wide  range  of  mesh  sizes.  Errors  over  1%  are  generally 

exhibited  when  the  ratio  of  the  partition  sizes  and 

on  C and  D,  is  greater  than  3 to  1.  The  error  sensitivity 

to  the  partition  ratio  of  H1/H2  is  understood  by  inspection 

of  Eq.  (3.2)  which  generates  the  influence  coefficients. 

The  dominant  coefficients  result  from  the  kernel  functions 

generated  within  and  between  the  C and  D crack  faces,  these 

1 2 

exhibit  a strong  dependency  on  the  values  of  (jj  ) as 

12  ^ 
compared  to  ( ~ ) . Only  off  diagonal  coefficients,  small 

H2 

compared  to  the  diagonals,  are  formed  from  the  kernel  func- 
tions interconnecting  the  circular  surface  and  the  crack. 
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In  general  these  problems,  and  the  bent  crack  geometries 
to  be  discussed  later,  require  the  restriction 
1/3  H^/H^  £ 3 to  insure  errors  of  less  than  one  percent. 

This  imposed  restriction  on  the  ratio  of  C/D  which  is  im- 
portant to  the  bent  crack  studies. 

Table  4.2  presents  example  problem  #1  results  for 
varying  partitions  and  increasing  crack  regions,  on  which 
the  unknown  displacement  functions  are  represented  by  the 
square  root  functional  form.  One  extreme,  assuming  linear 
displacement  behavior  over  the  complete  crack  boundary', 
shows  an  error  of  35%.  Rapid  reduction  in  error  (^2%)  is 
achieved  when  the  square  root  representation  is  utilized 
over  four  or  more  mesh  points.  These  results  confirm  the 
importance  of  properly  modeling  the  displacements  near  the 
crack  tip  and  provide  assurance  that  a reasonable  number 
of  mesh  points  leads  to  good  accuracy. 

The  data  in  Tables  4.1  and  4.2  illustrate  the  relair 
tive  insensitivity  of  the  results  upon  mesh  size  provided 
certain  guidelines  are  adhered  to.  This  results  from  pro- 
per modeling  of  the  unknown  displacements  and  careful 
numerical  integration  of  the  known  integrals  on  each  par- 
tition. It  is  most  important  to  model  the  square  root  be- 
havior of  the  displacements  in  the  neighborhood  of  the 
crack  tip.  The  accuracy  of  the  numerical  integration  scheme 
may  be  assured  by  successive  refinement  of  the  partitions 
or  use  of  higher  order  integration  techniques. 
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TABLE  4.1 


No.Eq. 

MT 

Ml 

M2 

M3 

Kl 

Cal. 

Kll 

Cal. 

Max.%  Error  in 
KI , Kji,  or  Es 

6 

5 

19 

.9979 

1.006 

.9983 

8 

7 

15 

1.009 

1.007 

.9964 

0.9 

98 

8 

7 

19 

1.009 

1.003 

.9991 

0.9 

106 

10 

7 

19 

1.009 

1.003 

.9989 

0.9 

114 

12 

7 

19 

1.010 

1.001 

.9989 

1.0 

130 

12 

7 

27 

1.002 

1.001 

1.0006 

0.2 

138 

10 

7 

35 

.9988 

1.002 

1.001 

0.2 

138 

12 

7 

31 

1.000 

1.001 

1.001 

0.1 

138 

16 

7 

23 

1.006 

1.001 

1.001 

0.6 

138 

10 

13 

23 

1.010 

.9991 

1.003 

1.0 

138 

14 

13 

15 

.9992 

1.006 

.9944 

0.6 

138 

12 

11 

23 

1.016 

.9993 

1.0022 

1.6 

138 

14 

9 

23 

1.013 

.9999 

1.001 

1.3 

138 

8 

15 

23 

1.024 

.9996 

1.003 

2.4 

Note:  1)  The  values  in  this  table  were  calculated  for 

C = .625,  D = .375,  E exact  = 1.561  #in/in. thickness 

s 


After  modeling  the  displacement  behavior  between  mesh 
points,  this  study  utilized  a Simpson  approximation 
coupled  to  an  accuracy  check  to  insure  errors  of  less 
than  1%.  Higher  order  integration  techniques  could  have 
been  selectively  used  to  save  modest  amounts  of  com- 
putational time  but  were  not  attempted  during  this  study. 
EXAMPLE  PROBLEM  # 2 . The  second  example  deals  with  the 
same  circular  geometry  except  that  a variable  length 
straight  crack  opened  by  a constant  internal  pressure  pQ 
is  considered.  J.  Tweed  and  D.P.  Rooke  (42)  report  re- 
sults for  a normalized  disk  of  radius  one  and  edge  crack 
length  a,  where  a is  allowed  to  vary  between  0 and  1. 

This  geometry  is  identical  to  that  shown  in  Fig.  7 with 

'Xj 

a = 0 and  C + D = a.  The  results  are  presented  in  Fig.  9 
for  values  of  a between  .025  and  1.9.  Comparison  with 
Tweed's  data  is  possible  in  the  range  ae(0,l)  and  shows 
less  than  a 1%  deviation.  The  results  of  Gross  et  al  (43) 
for  stress  intensity  factors  of  a single  edge  notch  speci- 
men, modified  to  represent  similar  internal  pressure 
loading  are  also  shown  for  comparison  purposes  in  the 
range  ae(l,1.8).  As  expected,  due  to  the  increased  stiff- 
ness of  the  straight  boundary  surface  vs.  the  curved  sur- 
face of  the  disk  geometry,  the  Gross  results  show  lower 
% 

Kj  factors  for  a values  of  0 to  1.4.  For  values  between 
1.4  and  1.8,  the  results  for  the  single  edge  notch  specimen 


65 


and  those  presented  here  for  the  disk  are  similar. 


The  results  of  Example  Problems  #1  and  #2  numerically 


confirm  the  solution  method  developed  in  Chapter  2 for 


plane  elastostatic  crack  problems.  In  addition  the  numerical 


techniques  necessary  to  obtain  accurate  solutions  are 


established. 


4.2  Bent  Crack  Geometry  Results 


The  last  example  considers  the  bent  crack  (a>  0) 


geometry  of  Fig.  7,  and  results  from  a parametric  investi- 


gation of  the  geometry  and  loading  are  presented.  The 


geometric  parameters  are  restricted  to  values  of  D = .5 


while  C e[.02,  . 3j  and  ae[o,85]  . A type  of  traction  load- 


ing is  examined  which  permits  a continuous  variation  from 


tension  to  shear  to  be  considered. 


Before  performing  the  bent  crack  study,  a straight 


crack  (ck=0)  with  C + D = .6  is  investigated  to  determine 


the  dependence  of  the  error  upon  the  ratio  C/D.  Results 


of  example  problem  #1  and  #2  (not  shown  explicitly)  are 


used  to  provide  a comparison  for  values  of  C/D  between  .01 


and  .3.  A partition  size  defined  by  Ml  = 10,  M2  = 19  and 


M3  = 27  is  found  to  result  in  the  lowest  error  for  the 


widest  range  of  C/D.  The  results  indicate  errors  of  less 


than  1%  for  C/D  >_  .05,  three  percent  error  for  C/D  = .025, 
and  eight  percent  for  C/D  = .01.  This  mesh  size  allows  a 


square  root  approximation  of  the  boundary  displacements 


on  six  points  in  the  neighborhood 


Mill  ipi* 


of  the  crack  tip.  As  discussed  in  section  3.2,  the 
piecewise  Y -displacement  behavior  is  incorporated  on 
10  mesh  points;  5 points  on  either  side  of  the  bend  vertex 
on  the  lower  crack  surfaces. 

A traction  loading  on  the  boundary  surfaces  which 
simulates  a continuous  variation  between  pure  tension  and 
pure  shear  on  a straight  cracked  disk  is  desirable.  This 
permits  comparison  of  data  with  reported  results  for  the 
bent  crack  in  an  unbounded  medium  subject  to  far  field 
loading.  This  simulation  is  accomplished  by  loading  the 
crack  faces  and  retaining  traction  free  boundary  conditions 
on  the  circular  domain.  A superposition  process  illus- 
trated in  Fig.  10  helps  to  relate  the  problem  with  trac- 
tion loading  on  the  crack  faces  to  one  in  which  the  cir- 
cular boundary  is  loaded. 

Begin  by  loading  the  circular  boundaries  of  an  un- 
cracked disk  with  tractions  F(P),  the  details  of  this  load- 
ing function  need  not  be  specified  at  this  point.  Assume 
that  on  the  line,  which  will  later  become  the  crack,  a 
constant  traction  of  magnitude  T and  orientation  P to  the 
x-  axis  results  from  the  applied  load  F(6).  Cutting  this 
line,  and  forming  a traction  free  crack  results  in  a crack 
opening  displacement.  The  crack  face  displacements  are 

brought  back  to  zero  by  the  application  of  crack  face 
_ 

tractions  - T.  Fig.  10a  shows  the  results  of  this  pro- 

. 

cess,  while  Fig.  10b  illustrates  a similar  cracked  body 
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with  crack  face  tractions  T applied.  Superposition  of 
the  two  problems  in  Fig.  10a ,b  result  in  the  problem  in- 
dicated by  Fig.  10c.  The  stress  intensity  calculated 
for  the  geometry  of  Fig.  10a  is  zero,  thus  the  stress  in- 
tensities for  the  geometry  and  loadings  of  Fig.  10b  and 
10c  are  identical.  The  strain  energy  release  rates  which 
may  be  calculated  from  the  stress  intensity  factors  will 
also  be  similar  for  the  geometries  of  Fig.  10b, c. 

Consider  a pure  Mode  I loading  by  defining  F (3)  to 
be  a constant  set  of  tractions  with  magnitude  T and 
orientation  3=  90°.  Constant  tractions  T applied  to 
the  faces  of  a straight  crack  (Fig.  10b)  leads  to  the 
desired  formulation.  Pure  Mode  II  loading  requires  F(0) 

IT 

to  be  taken  as  a set  of  tractions  Tr  = T cos  ( 2 6 — — ) 
and  Tq  = T Sin  (2  9+  where  9 is  measured  from  the 

x-axis.  This  results  in  a loading  for  the  problem  shown 
in  Fig.  10b  defined  by  constant  tractions  T applied 
parallel  (0=  0)  to  the  crack  faces,  and  opposite  in  direc- 
tion on  opposing  faces.  Constant  loading  on  the  crack 
faces  for  the  geometry  of  Fig.  (10b)  applied  at  angles 
0 < 3 <90°  result  in  mixed  mode  loading. 

The  geometry  and  loading  of  Fig.  10b  offer  further 
advantages  in  conducting  the  numerical  operations.  The 
higher  concentration  of  mesh  points  on  the  crack  surfaces 
compared  to  the  circular  domain,  coupled  with  the  simple 
loading,  produces  accurate  results  for  the  strain  energy. 


I 
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which  is  important  for  the  calculation  of  the  strain 
energy  release  rate.  As  a check  of  this  point,  a series 
of  straight  crack  geometries  under  Example  #2  loading 
are  solved  for  crack  lengths  of  .5,  .510,  .525,  and  .550. 
The  strain  energy  release  rate,  calculated  using  the  .5 
crack  length  as  the  standard  is  used  to  predict  the 
factor.  The  stress  intensity  factor  calculated  in  this 
manner  from  the  strain  energy  release  rate  is  accurate 
within  1%  when  compared  with  that  from  the  crack  face  dis- 
placement calculation.  This  calculation  requires  ac- 
curacy in  the  strain  energy  results  to  four  significant 
figures;  this  is  easily  achieved  for  all  of  the  check 
problems  considered. 

Results 

The  geometry  and  loading  of  Fig.  10b  is  investigated 
for  values  of  a between  0°  and  80°  in  5°  increments  and  $ 
between  0°  and  180°  in  15°  increments.  In  all  cases  con- 
sidered in  the  following,  the  magnitude  of  T is  set  equal 
to  unity.  Table  4.3  presents  the  calculated  values  of 
Ki'Kn  anc^  t^le  strain  energy  Eg  for  C = 0.25  and  D = .50. 
The  last  column  presents  results  for  a straight  crack 
length  C + D = .5  which  is  taken  as  the  initial  crack 
length.  The  expected  increase  or  decrease  in  Kj  and  Kjj. 
as  a function  of  a and  8 is  produced.  For  example,  at 
& = 0 which  represents  pure  shear  loading,  the  crack  at 

o = 0 begins  with  Kj  = 0 and  = 1.30.  As  the  crack 


TABLE  4.3 


BENT  CRACK  RESULTS 


'\a«= 

0° 

5° 

o 

O 

r—\ 

20° 

30° 

0 

.402879 
0 1.3 

.38513 
.25  1.02 

.371933 
.49  .82 

.372538 
.65  .72 

.373799 
.79  .59 

15 

.439216 
.613  1.25 

.431000 
.75  1.00 

.420990 
.91  .81 

.423354 
1.05  .68 

.426216 
1.16  .51 

30 

.538535 
1.18  1.12 

.530070 
1.22  .93 

.522079 
1.26  .74 

.525241 
1.37  .57 

.528815 
1.47  .40 

45 

.674223 
1.67  .92 

.66113 
1.59  .55 

.648113 
1.52  .63 

.650899 
1.61  .45 

.654104 
1.67  .26 

60 

.809923 
2.05  .65 

.78532 
1.89  .55 

.765322 
1.69  .47 

.766658 
1.73  .28 

.768513 
1.76  .10 

75 

.909274 
2.28  .33 

.872300 
2.00  .30 

.842301 
1.73  .27 

.841501 
1.74  .10 

.841386 
1.73  -.05 

90 

.945655 
2.36  0 

.898400 
1.92  0 

.858421 
1.67  .05 

.855373 
1.62  -.08 

.853196 
1.58  0.21 

105 

.909318 
2.28  -.33 

.86241 
1.83  -.08 

.809365 
1.48  -.15 

.804558 
1.4  -.26 

.800780 
1.32  -.35 

120 

.810000 
2.05  -.65 

.77108 
1.64  -.50 

.708276 
1.19  -.35 

.702671 
1.07  -.42 

.698181 
.98  -.48 

135 

.674312 
1.67  -.92 

.65234 
1.20  -.71 

.582242 
.83  -.54 

.577013 
.70  -.56 

.572892 
.56  -.56 

150 

.538612 
1.18  -1.12 

.50214 
.81  -.88 

.465032 
.41  -.68 

.461253 
.25  -.66 

.458483 
.02  -.61 

165 

.439261 
61  -1.25 

.40256 
.35  -.98 

' 

.388054 
-.04  -.78 

.386411 
-.20  -.70 

.385610 
-.34  -.62 

180 

.402879 
0 -1.30 

.35813 
-.25  -1.02 

.371933 
-.49  -.82 

.372538 
-.65  -.72 

.373799 
-.79  -.59 

Key:  E x 10  E = Strain  energy  lb-in/in. thick 

s s / 7 

Kj  KXI  Kx,  KII;~  psi  v in.  E = 10  psi 

D = .50  C_  = .025  v = .25 
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TABLE  4.3  (cont,) 


\«= 

p=\ 

0 

O 

50° 

O 

o 

70° 

**Q° 

0 

.375510 
.92  .43 

.377763 
1.04  .25 

.381188 
1.15  .02 

.408756 
1.32  -.28 

.362732 
0 1.1 

15 

.429538 
1.27  .34 

.433796 
1.38  .14 

.440637 
1.50  -.11 

.457678 
1.66  -.52 

.381135 
.56  1.1 

30 

.533054 
1.55  .21 

.538938 
1.64  0 

.559288 
1.72  -.27 

.576983 
1.90  -.71 

.467012 
1.1  1.0 

45 

.658319 
1.73  .06 

.665017 
1.78  0.14 

.678028 
1.84  -.40 

.724704 
2.0  -.88 

.584382 
1.5  .82 

60 

.771771 
1.78  -.07 

.778250 
1.8  -.26 

.792361 
1.84  -.51 

.833938 
1.97  -.97 

.701686 
1.9  .58 

75 

.843008 
1.71  -.21 

.848296 
1.69.  -.  37 

.861652 
1.71  -.59 

.902737 
1.81  -1.01 

.787603 
2.1  .3 

90 

.852944 
1.52  -.32 

.856387 
1.47  -.46 

.867334 
1.45  -.61 

.902666 
1.52  -.97 

.819071 
2.2  0 

105 

.793915 
1.23  -.44 

.800355 
1.16  -.50 

.807885 
1.11  -.61 

.833744 
1.13  -.86 

.787602 
2.1  .3 

120 

.695400 
.85  -.50 

.695213 
.76  -.52 

.699234 
.68  -.55 

.714439 
.66  -.72 

BB 

135 

.570134 
.42  -.55 

.569134 
.31  -.50 

.570494 
.21  -.46 

.576718 
.14  -.48 

.584450 
1.5  -.82 

150 

.456683 
-.04  -.55 

.455901 
-.17  -.46 

.456161 
-.28  -.35 

.457484 
-.38  -.24 

.467106 
1.1  -.10 

165 

.385445 
-.49  -.50 

. 385854 
-.61  -.36 

.386870 
-.74  -.21 

.388685 
-.88  0 

.381189 
.56  -1.1 

180 

.375510 
-.92  -.43 

.377763 
-1.04  -.25 

.381188 
-1.15  -.02 

.388756 
-1.32  -.28 

. 362732 

0 -1.1 

f 


bends  with  increasing  a the  Kj  factor  increases  as  is 

seen  to  decrease.  At  8 = 90°,  the  load  is  applied  normal 
to  the  straight  crack  and,  as  expected  at  o=  0,  K^.  = 2.36 
and  = 0.  Increasing  the  bend  angle  a results  in  de- 
creased values  of  Kj  with  commensurate  increases  in  K^. 
The  results  for  a constant  bend  angle  ot,  say  30° , as  a 
function  of  loading  angle  8 > are  also  of  interest.  As  the 
load  approaches  the  neighborhood  of  perpendicular  align- 
ment with  the  bent  crack,  for  example  8=  60°  for  a <=  30°, 
the  Kx  factor  and  the  strain  energy  release  rate  reach 
local  maxima.  As  the  load  angle  8 sweeps  past  this  neigh- 
borhood and  becomes  parallel  to  the  bent  crack,  crack 
closure  occurs  as  indicated  by  negative  KI  values.  For 
a = 30°,  this  first  occurs  between  8=  150°  and  165°.  The 
model  developed  in  this  study  includes  no  provision  for 
crack  surface  contact  under  compressive  loads;  thus,  the 
results  for  negative  K^.  are  of  little  interest. 

A careful  examination  of  the  strain  energy  data  in 
Table  4.3  for  8>120°  and  a > 5°  indicates  lower  calcu- 

lated values  than  the  corresponding  strain  energy  for  the 
straight  reference  crack  a=o.  This  indicates  a negative 
strain  energy  release  rate  contradicting  the  physics  of 
the  problem.  It  is  noted  that  in  all  cases  where  this 
contradiction  occurs  the  strain  energy  differences  are 
quite  small.  The  error  study  performed  for  ratios  of 
C/D  = .05  indicates  that  a 3%  error  might  be  expected  in 
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the  results  of  Table  4.3.  All  of  the  questionable  data 
in  Table  4.3  can  result  in  a zero  or  slightly  positive 
strain  energy  release  rate  when  a 2.5%  or  less  correction 
factor  is  applied.  The  results  would  then  be  consistent 
with  the  expected  physical  behavior. 

An  analysis  of  the  data  to  predict  the  initial 
direction  of  crack  growth  as  a function  of  load  is  readily 
performed  using  an  energy  criterion.  Given  a fixed  a for 
C = .025  and  D = .50,  the  resulting  strain  energy  as  a 
function  of  0 is  compared  to  the  results  for  the  straight 
crack  when  C + D = .5  for  the  same  3 load.  The  maximum 
strain  energy  release  rate  is  used  as  the  criterion  for 
determining  direction  of  crack  growth.  The  results  are 
shown  at  10°  increments  in  Fig.  11,  along  with  the  re- 
sults for  a bent  crack  geometry  in  an  infinite  medium  re- 
ported by  K.  Palaniswamy  (36)  . The  change  in  strain  ener- 
gy as  a function  of  a and  3 is  found  to  have  a relatively 
flat  maximum  extending  over  a range  of  10  to  15°.  The 

sharpest  maxima  occur  for  values  of  a near  zero  and  3 near 
o 

90  . The  data  in  Fig. 11  have  bars  indicating  an  estimate 
of  the  error  incorporated  in  numerically  selecting  the 
critical  values.  These  flat  maxima  have  been  studied  by 
J.  Swedlow  (45)  for  similar  problems  of  interest.  He 
appropriately  concluded  that  small  changes  in  local  micro- 
structure could  strongly  influence  the  growth  direction 
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making  experimental  study  difficult.  He  presents  both 
analytical  and  experimental  evidence  to  confirm  this 
behavior. 

The  numerical  results  for  the  bent  crack  problem  in- 
dicate the  presence  of  a critical  angle  ot^.  This  critical 
angle  manifests  itself  through  the  geometric  influence  co- 
efficient matrix  [C]  which  develops  a zero  determinant  for 

a = a and  results  in  large  values  of  K_,KTT  and  E as  a 
c ^ I'  II  s c 

is  approached.  The  geometry  reported  in  Table  4.3,  rep- 
resenting a C/D  ratio  of  .05,  shows  that  aQ  equals  78° 

(to  the  nearest  degree) . Lower  values  of  C/D  lead  to 
similar  results  while  larger  C/D  values  produce  lower 
critical  angles.  These  results  are  summarized  below. 


c 

D 

C/D 

a 

.015 

.50 

.03 

78 

.020 

.50 

.04 

78 

.025 

.50 

.05 

78 

.050 

.50 

.10 

63 

.100 

.50 

.20 

58 

.150 

.50 

.30 

46 

.200 

.50 

.40 

42 

.300 

.50 

.60 

38 

.400 

.50 

.80 

35 

It  should  be  noted  that  the  numerical  value  of  the  deter- 
minant for  all  problems  is  quite  large  (10^  to  105*"*)  . 
This  determinant  remains  relatively  constant  for  values 
of  a<  (ac  - 5°) , then  decreases  rapidly.  Once  past  the 


I 

■* 


critical  angle  otc  the  determinant  assumes  large  negative 
values.  The  calculated  for  a>  ac  becomes  negative, 
thus  the  results  for  K_,KTT  and  E , for  a>  a are  invalid. 

' As  background  to  this  problem  recall  the  work  of 
Sternberg  and  Koiter  (39)  . Their  studies  for  an  open 
wedge  geometry  subjected  to  a concentrated  moment  at  the 
tip  focuses  on  the  break  down  of  the  solution  for  wedge 
angles  of  257.45  degrees  (this  corresponds  to  a = 77.45° 
for  the  bent  crack  geometry).  Westmann  (40)  showed  that, 
at  a=  77.45°,  a second  singular  contribution  may  enter 
the  analysis  while  only  one  singularity  exists  for  a>77.45. 
His  work  centers  on  the  geometric  effects  in  adhesive 
joints  and  illustrates  the  importance  of  this  second  sin- 
gularity for  both  the  development  and  elimination  of  stress 
intensity  factors  at  the  wedge  vertex.  K.  Palaniswamy  (36) 
discusses  the  existence  of  a critical  angle  for  Mode  II 
loading  and  reports  a value  of  77.4°.  He  does  not  enter 
into  detailed  discussion  but  indicates  the  existence  of 
qualitative  experimental  results  taken  in  the  neighbor- 
hood of  a . 

c 

Certainly,  the  existence  of  ac  as  at  least  a mathe- 
matical phenomenon  is  well  documented.  However,  the  work 
of  Westmann,  Palaniswamy  and  experimental  observations 
for  Mode  II  loaded  cracks  leave  open  to  question  whether 
a c is  indeed  a physically  real  limit  upon  the  crack  growth 
direction. 
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The  resulting  lower  values  indicated  for  as  the 

bent  crack  C extends,  are  in  contrast  with  Palaniswamy ' s 
results.  Since  he  reports  in  Fig.  12  of  (36)  increasingly 
permissible  bend  angles  for  increased  values  of  C,  the 
bent  crack  length.  The  results  presented  in  (35)  are 
obtained  for  a bent  crack  in  an  unbounded  domain  while 
the  problem  considered  here  is  for  a bent  edge  crack  in 
a finite  disk.  The  contrasting  differences  in  the  two 
results  may  be  caused  by  crack  tip  interaction  with  the 
boundary  and  severe  bending  induced  for  larger  C/D  ratios, 
features  not  present  in  the  problem  considered  by 
Palaniswamy. 
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CHAPTER  5 


CONCLUSIONS  AND  FUTURE  APPLICATIONS 

The  results  presented  herein  demonstrate  the 
application  of  integral  equation  analysis  to  two-dimen- 
sional elastostatic  fracture  mechanics  problems.  A solu- 
tion method  is  developed  for  analysis  of  general  closed 
edge  cracked  bodies  subject  to  traction  loading.  This 
method  is  demonstrated  using  example  problems  with  known 

solutions,  the  results  for  KT , K„._,  and  E beinq  accurate 

J.  II  s 

within  one  percent.  The  bent  crack  geometry  illustrates 
the  application  of  the  method  to  an  asymmetric  complex  do- 
main which  incorporates  the  effects  of  singularities  at 
the  bend  vertex. 

The  application  of  this  solution  method  to  more  gen- 
eral two-dimensional  edge  cracked  problems  is  possible. 
However,  the  practical  limitations  of  the  method  should  be 
clarified.  Section  2.3  discusses  the  solution  for  a gen- 
eral plane  cracxed  body.  In  the  closing  portion  of  Sec- 
tion 2.3  it  was  implied  that  a body,  with  m external 
cracks  and  N parametric  relationships  describing  them,  will 
require  a (4N  + 2)  system  of  coupled  integral  equations  to 
be  analyzed  and  solved.  Therefore,  the  bent  crack  requires 
ten  such  equations.  Geometries  with  non-circular  outer 
boundaries  requiring  multiple  parametric  representation 
(say  N2)  lead  to  further  expansion  of  the  equation  set. 


The  formula  determining  the  size  of  the  equation  set  now 
becomes  4N  + 2^.  » For  example,  the  simple  rectangular 
plate  with  one  bent  edge  crack  requires  five  parametric 
expressions  for  the  outer  boundary  and  two  for  the  crack 
faces;  a system  of  18  coupled  equations  results.  Thus, 
analysis  of  complex  geometries  by  the  method  represented 
herein  may  require  manual  manipulation  exceeding  the  prac- 
tical value  of  the  solution.  Future  developments  in  com- 
puter software  technology  may  alleviate  this  difficulty. 

It  is  currently  possible  to  construct  programs  such  that 
the  geometric  coefficients  resulting  from  most  of  the  equa- 
tions are  automatically  calculated  from  simple  boundary 
surface  data.  The  computer  derivation  of  the  equations  re- 
sulting from  the  limiting  process,  -►O  and  differentia- 
tion, w.r.t.  ij>  is  beyond  the  current  computer  software 
technology.  Successful  development  of  this  technology  will 
permit  practical  application  of  this  method  to  complex 
geometries  without  extensive  problem  preparation. 

For  the  present,  analysis  of  complex  structures  us- 
ing finite  elements  is  possible.  This  method  incorporates 
special  singular  elements  adjacent  to  the  crack  tip  and 
provides  the  analyst  with  an  accurate  tool  for  the  deter- 
mination of  stress  intensity  factors  in  complex  geometries 
(17,46-50).  The  circular  element  with  bent  crack  developed 
in  this  study  may  possibly  be  incorporated  as  the  special 


element  enclosing  the  crack  tip.  This  will  allow  study  of 
complex  geometries  with  bent  cracks, and  at  the  same  time 
incorporate  all  of  the  advantages  of  the  finite  element 
method . 

The  success  achieved  using  integral  equations  in 
the  analysis  of  two-dimensional  elastostatic  problems  led 
Cruse  (21,22,51)  to  the  natural  extension  of  the  solution 
capability  to  three-dimensional  problems.  In  (51)  he  pre- 
sents the  B.I.E.  for  three-dimensional  problems  valid  for 
reference  points  P(x)  not  located  at  an  edge  or  corner. 

His  equation  (14)  yielded: 


i u.(P)  + ui(Q)  T..(Q,P)  dS 
S 


ti(Q)Uij(Q'P>dS(Q) 


(5.1) 


where : 


0,,  . i—  i 6..,  + 1 


'ij  4 iry  r 4(1  -v ) ’ 4(1  -v)  r'ir'j 


K = (1  - 2v)/2(l  -v) 


Tij  -4tt  r2  3 n (<Sij  + 1 -2V  r'ir,j)  " njr'j+njr'j 


where  the  integrals  are  interpreted  in  the  sense  of  the 
Cauchy  principal  value  and  the  notation  is  similar  to  that 
in  Chapter  2.  Eq.  (5.1)  can  be  viewed  as  the  governing 
equation  relating  surface  tractions  to  surface  displace- 
ments. Eqs . (2.14)  present  the  equations  for  the  two- 


dimensional  problem.  The  mathematical  character  of  Eqs . 
(2.14)  and  (5.1)  is  similar,  leading  to  the  speculation 
that  the  solution  method  developed  and  verified  for  anti- 
plane and  plane  two-dimensional  problems  can  be  extended 
to  three  dimensions. 

In  summary,  this  work  extends  the  B.I.E.  developed 
by  Rizzo  to  the  application  of  closed  crack  geometries. 

It  demonstrates  that  accurate  numerical  solutions  can  be 
obtained  by  this  method  for  finite,  asymmetrical  cracked 
bodies . 


* 
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APPENDIX  A 


GEOMETRIC  KERNEL  FUNCTIONS 

This  appendix  presents  the  geometric  kernel  func- 
tions for  Eq.  set  (3.1)  which  apply  to  the  geometry  of 
Fig.  7.  The  one  hundred  kernels  are  presented  in  twenty- 
five  sets  of  four.  Each  set  has  a heading  Ker  L M where 
L,  M have  values  1 tc  5.  The  L index  defines  the  reference 

surface  PT  and  the  M index  defines  the  field  surface  I\_ 

L M 

to  which  the  kernel  corresponds. 
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Tsin0-  C- tana+B’ sin (a-^n)l 

= K-M’4*vr 


• ri-E’cos0-C’tr-na’Sin0-B’cos(0+a-<(»on/R2  1 

ll  "H/ 


Ker  21 


„ f 2 2 ] 

R2  = <,  (E+p  •cos(a-i/;o)-B2)  + (C‘  tana-p  • sin  (a-iJ>o) -C*  tan’J'  ) ^ 

~21  fr  2i  ° ‘ 

Kl2i=  K-M*4*p  -~E— P- co s ^ a to)..,5?)  [(E+p  • cos  (ct-^o)-B2>  (-sin  (ct-tpQ) ) 

L R2  J 

v « i 


+ (C* tana-p -sin (a-^n) “C'tan^o) (-cos (a-^n) )1/R2  ; 

“21 J 


K1I2i=  |“M*4*p[e  +p -cos (a-^g)  -b2]§:‘ tana-p  • sin (a-^g) -C*tani|»o]  • 

[[(E+p  • cos  (a-i|ig)  -B2)  (-sin  (a— v/jq  ) ) 

+ (C*  tana-p  • sin  (a-i^g)  -c*  tan^o)  (-cos  (a-iJiQ ) )]  /R^ 

+ ^2E<E+P<:os(a“'l'o)-B2Hcos(a-((;o)i 


+ (C*  tana-p-sin  (o-iJiq)  -C’tanifiQ)  (-sin  (a-iJiQ) )] 


KHI2i=  |^-M* 4* p [e+p  - cos (a-^o ) -B2]  [c- tana-p  • sin  (a-^0)  “C'tanijio]  ' 

[[(E+p  ’cos(a-i^o)  -B2)  (-sin(a-ii<o) ) 

+ (C-tana-p  •sin(a-i/jn)-C'tan4io)  (-cos(a-i|»o) Jj/R1* 

~2 1 

- p[(E+(Kffls(a-i|)o)-B2)  (cos(a-ipo)) 

^ 1 

+ (C*  tana-p-sin  (o-^q)  -C'  tamj/p)  (-sin(a-^o) )]  i 


kiv2i=  s[k-m*4*u 


(C* tana-p  ’sinfo-iZ/nJ-C-taniJ 


[(E+p -cos (a-tyo ) -B2 ) (-sin(a-^o) ) 


+ (C-tana-p  • sin(a-t)io) -C-tan^n)  (-cos  (a-^p)  )]/R2  ! 

“2 1 1 


91 


Ker  23 


2 2 
R2  = (cos0-B2)  +(sin0-C'tanij)o) 
~ 23 


KI23=  K - 


M*  4-  p (COS0-B?) 


> 

-jjVP-Zl 


cos0-  C'tamjin  • sin0 


KIl23=  l -M*4"p  [cos0-B2]  [sin0-Ctaiu|>oJ 


□L-Bo  * cos0-  C*taml>n  ’ sin0 


+ K- fB2*sin0-  C- tcini|)o -cosOl/R2 

“23 


r 

KIIl23=  j^-M*4* p [cos0-B2]  [sin0-C'tan\j>g] 


cos0-  C-tanjJ 


K*  [b2*  sin0-  C-tani|>o'cos0j/R^ 


KII24  = ‘ 

j-M*4 

_ J 

,-p[xi-b2 

+(xi-b?) 

] [-2*C*  tan^o]  2/R4 
~ 24 

]) 

"1 

~24 

KIII24  = * 

i 

l-M-4 

L 

+ K 

,*u  [xT-b2 
xT-b?] 

R2 

L "24  J 

] [-2*C*tan4io}  /R4^ 

> 

+ K 


R2  = (E+p*cos(a+^g)-B2)  + (C* tana-p -sin (a  +tj>g)  -C*  tan>|>g) 

~ 2 5 

l 

ki25-  |[k-  w~-4'-"  1e+P  [lEtp.eo.(eri.<,|,,-B2)  si„«,+*0> 

R2 

l “25  1 

. I 

+ (C*  tana-p*  sin  (a+iJig)-C'tan^g)  cos (a+ij>g) ] /R^ 


KIl25=  ^-M-4*p[E+p*cos(a+ij;g)-B2]  [c*  tana-p  * sin  (a+ij>g ) -C*  tantjjg] 

■ [(E+p*cos(a+i|(g)-B2)  sin(a+ij;g) 

+ C*  tana-p  • sin  (a+ijig)  -C*  tani/ig][cos  (ot+^g ) J /r£ 

+ ~2  C<E+P  * cos  (a+ij^g ) -B2 ) (-cos  (a+ijig) ) 

+ (C*  tana-p  * sin  (a+i|>g ) -C * tantjjg ) sin  (a+ijig ) J 


KIII25—  <j~M*4* p [e+p *cos (a+$g) -bJ  (c* tana-p  • sin (a+tpg ) -c*  tanipgj 
• [(E+p-cos(a+i|»g)-B2)  sin (a*H|»g) 

+ C* tana-p  • sin  (a+i|>g) ~C*tani|;g]  [cos(a+tyg)]/R^ 


■ ^2  [(E+p-cos(a+i|)g)-B2)  (-cos(a+<»g)) 

+ (C* tana-p*  sin(a+t|)g)-C'tant|)g)  sin(a+tj>g)] 


lr.  M*4*p  (C* 

KIVo  5 = \ K-  

a 

[_(E+p  • cos  (a+i^g ) -B2 ) sin  (a+i^g ) 


tana-p • sin (a+fln ) -C* tan^n ) 


“25 


+ (C*  tana-p  • sin(a+^g)  -C*  tani|>g)cos  (a+tjig)]]  /R2 

~ 2 5 


95 


Ker  31 


1 


'31 


(E+p*cos(a-t(/o)-cos0o)  + (C*tana-p*sin(a-iiio)-sin0o) 


KIV31 


1 


K-M*  4*  p (C*  tana-p ; sin (a-»n) -sin0n ) 


"31 


[_-E*  sin(a-\|»o)-C*tema'cos(a-i|/o)+sin  (0o+a-\Jio)]/R^ 


ki3i= 

« 

•[- 

2 

fir  m a-,,  (E+p*cos(a-(pn)~cos0n)  -1 

1 

Ln(0o+a-i|)O)]/R2  \ 
“31 

J 

L r2 

. 11  J 

•E*  sin  (a-i)»o)  **C*  tana*  cos  (a-<J»o)  +sj 

KII3i=  < 
+ K*  [e*  COE 

r 

-M* 4 * p [e+p * cos  (a-t|)o)-cos0ol  [c*  tana-p  • sin  (a-iJiQ)  -sinOp] 

:*sin(a-i|)0)-C*  tana*  cos  (a-\J»o)  +sin  (0n+a— 1^0)]]  /R1* 

31 

; (a-^o)  -C*  tana*  sin  (a— t/»o ) +p-cos  (0Q+a-^o)]  /R2  | 

KIII3  j=  ^-M*  4*  p [e+P  * cos  (a-1^0 ) “Cos0o]|  [c*  tana-p  • sin  (a-ijip) -sin0o] 

• [-E*  sin  (a-\|»o ) -C*  tana*  cos  (a-i|)Q ) +sin  (0Q+a-iJi0  /R4 

31 

- K*  [E*cos(a-\|io)-C*tana*sin(a-«Jjo)+P**cos (0o+a-iJio)J/R2  l 

~31J 

kii32  = •- 


+M*4*u[xi-cos0ol  f C'tani|/o-sin©ol  /R4 


+ K • [XI  -cosGg]  /R2 


'32 


[ 

KIII32  = ) M-4-y  (xI-cosGo]  [ C-tani{;o-sin0o]  /R4 


r 


32 


K • [xi-cosG0]/R2  > 

~32  I 


Ker  33 


* 


R2  = 2 • ( 1-cos (0-0o )) 

~33 


KI33  = < 
•( 

1 

|j-K_M.4.pJcos0z_cosMij  j-  i j 

■ 

► “P  [■*»<•••.>] + 1 • =fs£-! 

sin(  ) 

! 

! M-4-y  K cosf^11) 

|+  4 -Sin(0+0O)  - 2 • 

; 

kxv33= 

h -MT  • cos^  i6;0”)} 
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Ker  34 


R2  = (Xl-cos0o)  2 + (-C*  tant|/Q-sin0o)  2 


( 

ki34  =| 

r m 1 (Xl-cos0n) 2 

-sin0n-C*  taru|>r> 

' 

R2 

“34 

R2 

34 

KII34  = . 

f 2 

j-M*  4-  p ]xl-cos0o]  [C-  taruJ/o+sinOgl/R4 

l | 34 

+ K ■ [cos0o-Xlj/R2  > 

34  j 

KIII34  = - 

2 

-M*4-p(xr-cos0Q]  [c- taruj>o+sin0o]  / R4 

- K • [cos0o-xTj  /R2  j> 

^ .J 

f 

r .r- 

KIV 


34  ~ 


K-M'  4 • p 


(C~  tarn|)n+sin0n)  2 

R2 

34 


-C,taml)n~Sin0n 

R2 

34 
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Ker  35 

R2  = (E+p‘cos(a+iJig)-cos0g)2  + (C*  tana-p  • sin  (a+ijig) -sinOg) 2 
35 


KI35  = JfK-M-4-p  P - COS  (ct+^0 ) ~,cos0fl ) : 


L35  = 1 

R2 

"35  J 


* 


[jE*  sin  (a+<|>g)+C*  tanacos  (a+i|ig)  -sin  (a+\Jig+0g)]  /R2 


KII35  = ^ -M-4-p[E+p-cos(a+i|ig)-cos0g]  [c ‘tana-p*  sin  (a+i|ig) -sin0gj 

• [E*sin(a+iJ;g)+C*tanacos(a+iJig)-sin(a+\Jig+0g)]/R^ 

+k[j-E  • cos  (a+i|;g ) +C*  sin  (a+i|ig ) • tana-p +co  s (Og+a+^g )]]  /R2 


35 


KHI35  = | -M-4-p[E+p*cos(a+<Jio,  cos0g]  [c*  tana-p  • sin  (a+\Jig)  -sin0g] 

• [e*  sin  (a+<|;g ) +C’  tanacos  (a+<|>g ) -sin  (a+ifig+0g)]  /R4  , 

3 5 


-k[-E-  cos  (a+<|;g)  +C‘  sin  (a+i|/g ) * tana-p+cos  (0g+a+\Jig)]  /r‘ 


35 


Ker  41 


R2  = (E+p  *cos  (a-^o)  "82)  2+  (C.  tana-p  • sin  (a-\J;o)+C*  tan'i'o) 2 
~ 41 


f 


Kim  = 


K-M-4-y  SEtBl£2S^ial^2)  : 


"1+1 


[(E+p  * cos  (a-^ig ) “B2 ) (-sin(a-'i'o) ) 


1 


KII41  = 


+ (C*  tana-p  * sin (a-^g) +C* tan^o ) (-cos (a-i^g) )]  [/r2 

! 4 1 


-M*4-p  [s+p  *cos  (a-^o)  -82]  Cc‘  tana-p  • sin  (a-i^0)  +C*  tani^oj  1 


Qe+p ’ cos (a-ij)g ) -B2 ) (-sin(a-^0J) 

+ (C* tana-p  -sin  (a-iJjrO+C*  tan^g)  (-cos(a-ij'g)  )J /r4 

— k 1 

K 

R2 
~kl 


+ -2[(E+p-c°s(a-i|)o)-B2)  (cos  (a-i|>g )) 


+ (C*  tana-p  • sin  (a-ij>g ) +c‘  tanijig ) (~sin  (a— 4»0 ) )] 


KIII41  = j -M* 4* p[e+p- cos (a-iJigJ-Bj)  tana-p* sin (a-i(ig)+C*tan^g] 

[jE+p  • cos  (a-ipg ) -B2 ) (-sin(a-ipo) ) 

+ (C-tana-p  -sinta-ifigJ+C’taniiJg)  (-cos  (a— ) )]/r4 

h] 

" §2  [(E+p*cos(a-i|>0)  -b2)  (cos(a-ijjg) ) 

_1+1  1 

+ (C*  tana-p  • sin(a-ii>g)+C*tani()g)  (-sin(a-ipg)  )J 


KIV4 


1 = |[k  - M*  4 ’ p — 


tana-p  • sin  (a-ipn ) +C*  tan^n)  ‘ 


~41 


Q(E+p  • cos  (a-t|)0 ) -b2  ) (-sin  (a-i//0 ) ) 

+ (C*  tana-p  • sin  (a-\J>g)+C*  tani^g)  (~cos  (a-i|>g) 
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U 


Ker  42 


R2  = <(X1-B2)z  + (2*C-tan\|i0)‘ 


“42 


KI42  = 


K - 


M-4-y (Xl-B?)‘ 


~ 42 


-2‘C«tan^n 


“42 


kii42  = 


+M-4-W (Xl-B2)  (2-C*tamJ<o)2/R4 
' “42 


+ K 


XI  - B? 


R2 
-42  J 


KIII42  = \ +M* 4 * p (Xl-B2)  (2-C*tani|)0)2/R4 

| ~42 


- K 


XI  ~ B? 
R2 

"42  J 


KIV42  = 


T,  rf4-u(2*C‘tan<'n)2 

-2*C*  tan<J>n 

' 

K-M  4 M 

R2 

R2 

“42 

“42 

~ 

j 

102 


= (cos0-B2)2+  (sin0+C- taring) 2 


KI 


43 


[k-M-4-u  ^cosQ~b,>)2  I rl~B?-  cosQ+C-  tanifrn*  sin0 


“43 


- rl 


“43 


= j-M-4-1 


KII43  = |-M-4-u(cos©-B2) (sin0+C-tan*o) 


(l-B? •cos0+C‘ tan^n  * sin0) 


+ K[B2*sin0+  C- taring* cos0l  /R2 

“43 


"43 


KIII43  = < 


-M*4-p(cos0-  B2)  (sin0+C‘ tanipg) 


K [b2 - sin0+  C- tamin' co s0]  /R2  , 

"4  3 


(1-B?*cos0+  C'tan^n* sin0) 

R4 

"43 


KIV43  = 


f7r  M * „(sin0+  C-tan^n)21 

r 1-B9‘cos0+  C- tanij/n  • sin0  -|  | 

L R2  J 

( 

L "43  J 

L “43  J J 
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. i-.-—  .-.',.*1*- 


i 


Ker  44 

R2  = (XI  - Bn)2 
“44 

KI44  - 0 

Kill,  4 = — 

XI  - B2 

KIII44  = ~ — 

XI  - B2 

KIV44  = 0 


104 


Ker  45 


R2^  = (E+p  • cos  (a+i|>o ) -B2 ) 2 + (C- tana-p- sin  (a+i|>g)+C*tani|>g)  2 


= [k-M-4-u— +P-‘- 


1 _n„)2 


hr  • cos  (a+i|>g ) -B2 ) sin  (a+\|>g ) 


+ (C- tana-p  • sin  (a+\|>g ) +C*  tani|>g ) (cos  (a+\|>g ) )]  /R2 

~ U5 


KII45  = s -M- 4- y (E+p • cos (a+i|>g ) -B2)  (C- tana-p • sin (a+t|>g ) +C-  tan\|>g ) /R** 

l ~**S 

• [(E+p  ■ cos  (a+\J>g ) — B2)  sin(a+t|>g; 

+ (C*  tana-p  • sin  (a+\|>g ) +C-  tan^g ) (cos  (a+\|>g ))] 


+ ^2[(E+p-cos(a+\|>g)-B2)  (-cos(a+^g)) 


+ (C  • tana-p  • sin  (a+\J>0 ) +C  • tan^g ) sin  (a+i|>g  )J 


KIII45  = M*  4 * u (E+p  • cos  (a+i|fg)  -B2)  (C*  tana-p  • sin(a+i|>g)+C-tani|>g)/Fr^ 

• [(E+p  - cos  (a+\|>g ) -B2 ) sin  (a+\|>g ) 

+ (C-  tana-p  • sin  (a+\|>g ) +C  * tani|>  g ) cos  (a+i|*g ) ] 

- |2  [e+p  •COS(a+^g)-B2)  (“COS  (a+l|>g)J 
~4  5 

+ (C- tana-p  • sin  (a+\|>g)+C-tan\|>g)  sin(a+\|>g)J 


KIV45 


=1  Ik-  m-4-p— 


tana-p • sin (an 


| (E+p-cos(a+\|>g)-B2)sin(a+\|>g) 


+ (C-  tana-p  • sin  (a+t|<g ) +C-  tantjig)  •cos(a+t|>g)]/R2 

4 5 
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* 


Ker  51 


R2  = p2+B2  -2*pB*cos  2ipn 
“51 


KI51  = 


. (p*cos{a-^n)-B*cos<a+i|>n))2i  r 

' 

\i  j ( 

R2  JL 

~51  L 

** 

-5.  Jj 

[p*cos(a-ifio)  -B-cos(a+ipo)] 

B*  sin(a+ipo)]  . T-?— + K 

rP-B-cos  2i|>n 

L R4 

L R2 

I 

~5l 

L ~51 

4 

KIIS1  = 


KIII51  = s 


M*  4*p  [p "cos (a-ipg)  ~ B-cos(a+^o)] 


l | 

• [>-.ln(a-»i,)-B-«ln<Mft)3  • rB'sl"  ^"l  - xf  ‘>~B'cos  *♦»  1 

L 4 j L 4 Jj 


kiv51 


K-M*  4*  p 


(p-sin(a-*|Jn)~B.sin(ci+il>n)  )2  t r -B*sin  2^0  -i  j 

R2  JL  R2  J 

C 1 -J  I—  — f 1 -j 
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»»»»>  winii  n-i'niM'drX  fc  » ■ .» . 


Ker  52 


l]l  UIJI  u II  , I «T 


— = [xT-E-B'Cos  (a-Hf/g)]  2 + [c*  (t«mi(io-tana)+B-sin(a+^o)J  2 


ki52  = | 

r..  M-4-u  (xl-E-B'cos  (a+ii»n) ) 2 T 

“52 

1 

• [^-C(tani(/0-tana)-B- sin(a+\(;o)J/R2  f 

KII52  = < 

: 

fxi-E-B.cos(a+^)l  * [C-  (taniJ;n-tana)+B-sin(o+iJ(n)l2 

R4 

52 

+ K * [>a-E-B-cos(ci+^0)]/R  2 > 

52 

kiii52  = 1 

M-  - p fxi“E-B- cos (a+ipn )1  * fc*  (tani{»n-tana)+B-sin(a+i|>n)l  2 

l a4 

521 

- K • [xi-E-B  • cos  (a+ipo ) J /R2  ? 

■"J 

f 

kiv52  = . 

• 

rK  h-4...  (C- (tan\|jD-tana)+B-sin(a+iJ)n))2  7 

L R2  J 

-52 

C-C.  (tanip0-tcina)-B’sin(a+i|)o)]/R2 

S2  f 

I 
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H. 


Ker  54 


R2  =^Xl-E-B*cos(a+i|/o)]2  + Q-C‘taniii()-c‘tanot+B'  sin(a+i)Jo)J  2 
54 


M‘4*yi  (Xl-E-B* cos (ct+ijin  H)2-!  r~C‘  tanifrn-C-  tanot+B# sin ) 


KI54  = K- 


“IP 


KII54  = j -M*  4 tj|  Xl-E-B*  cos  (a+iI)n)1T-C‘  tan^n-C'tana+B*  sin(a-n^n) 


- K • [xi-E-B-cos(a+i|>0)]/R2^ 


KIIj5  _ I M ^.^[xT-E-B-costa+iiJn)!  • r~C<tan^n~C-tana+B-sin(a+il)n) 


+ K • [xi-E-B-cos(a+iJj0)]/R2  , 

54 


kiv54  = 1 fie  - M,4‘Ur-C-tani(in-C-tana+B-sin(a+0n)l2l 

I R2  J* 

L L “54 

|-  -C'tani^n-t 


C ' tanij;  n -tana+B  • sin  ( a+i|)  n ) - 


F 


APPENDIX  B 

GOVERNING  B.I.E.  FOR  CLOSED  BENT  CRACK 

This  appendix  contains  Eq.  set  (3.2),  the  governing  B.I.E.  for  the 
closed  bent  crack  (^Q  * 0).  Section  B.l  presents  the  full  set  of  equa- 
tions. Section  B.2  presents  algebraic  formulas  inserted  into  these 
equations  depending  on  the  assumed  displacement  model.  Section  B.3 
contains  the  detailed  forcing  function  formula  necessary  for  general 
traction  loading. 


Define : 


R2  =|xl  - E - B • cosa]2+  [-C  • tana  + B • sinaj- 
12 

R2  = Tsin0  - E - B - cosaj  2+  [sin0  - C • tana  + B 
13 

r2  =|e  + p * cosa  - B2]2+  [c  • tana  -p  • sina]2 
21  u 

r2  =0.-2  • B2  • cos©  + B2~] 

23  l" 

R2  = [e  + p • cosa  - cos0gl 2 + [c  • tana  - p • si 
31 

r2  = fl  - 2 * Xl  • cos0q  + Xl  2J 


x [U1(B)J  + I U2  (XI) 


M • 4 • P (Xl-B-B- cosa) 


, 2ir 

J U3(0)  [K-M‘4' 

0 lL 


u (cosQ-E-B^COSg) 2 1 

R2  J 

13 


p.-E*cos0-C‘tana'sin0-  B* cos (0+q)] /R2  d0 

r1  I r M-4-urxi-E-B-cosa] 2 C-tanq-B-sinq' 

' I U4  (XI)  j R2  R2 

-p  I L 12  H 12 

,B-H 

’ I • t>5(B)]  + J Vl(p)‘  ^b”  dP 


+V1(B-H)  *(A1VMH)  +V1(B)  ’ U1VB)  +V1(B+H)  • (A1VPH) 


Vll 


♦ vl(pi'T:r  dp 


) . fxl-E-B’cosql  r-C-tanq+B- sinql 2 

+ v — 1M*4'iil J 

2 (XI)  | 11  R4 


+ K [Xl-E-B*  COSqJ  /R2  | <3X1 


3.2.A.1 


3(0) 


| -M*4*y  [cos0-E-B-cosa],[sin0-C*teina+B‘  sina] 


fl-E-cosO-C’ tana  * sin0-B‘cos  (0+a)]  /R1* 
L 13 


K • [E*sin0-C* tana ’cos0+B‘ sin (0+a) J/Rj3d0 


4 (XI) 


, , [xi-E-B'cosal-[+B- sina-C‘ tanal ' 

+M*4*y  i-= J 


12 


+ K [xi-E-B*  cosa]  /R2 


12 


B-H 


dXl 


5(p)  p-B 


dp 


• V5(B-H)  • (UVMB>  - V5(B)  ■ (ilVBl  - V (B+H)  •(“VPB> 

(£ 

'J  v5(p> * w dp  " I ' FM(B) 

B+H 
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hJ  V2(i)  {E?i2}  dX1  +j  v3(0)  |-m-4'u[cos0-B2] 


B2+H2 


&ine](i^sai]  . K . [B2.sin9]/R23) 


*\  2 ‘xTS-,1  dX1  - - V4(B)-<S2VB> 


-v4(b,  MAOTH,  +j  V4(-  dXl 

b2+h2 


+ j v5(p)  [e+p  •cosoi-b2]-[c*  tana-p  • sinaj 


Re+p  • cosa-B2)  sina  + (c* tana-p  • sina) co sal /R4 

21 


+ K Re+p  • cosa-Bo)  (-cosa)  + (C- tana-p • sina) sina] /R^  Idp  = 7-  • FA2  , . 

21|  i.  (B2) 
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1 

I 


117 


dXl 


[\-J 

J 2 (XI)  < 


M.tj;5r-cose0]fei2?J2+  k • taLcosjj 


32 


32 


•/'■  - 


3(0) 


M“4"U  • 

— -sinO+Qo)  j d@ 


0o-h3 


K cos  (~2~^") 


'3(0)  1 2 7”B-U 


sin 


d0 


+ K • { 

l 


V3(0q+H3)  V3(0q-H3) 


2it 
+ I v 


j 


3(0)  2 


0O+H3 

a 


K [ cos  (^r^4 


4 (XI) 


r—  „ l&inQnl2  „ [xi-cos0nl 
-M"4-y j_xl-cos0ojc;t^yu-  - K*1 ^7 — 

32  3? 


> dXl 


( 


t 

V„ . * j -M-4* 


5 (p) 


p[e+p  •cosa-cos0o]*[c-tana-p • sina-sin0o] 


[e- sino+C* sino-sin (a+0o )] /R^ 


+ K • £-E-cosa+C’sina’tana-p+cos(0o+co3/R^  l dp 


— • FA3 

X FAJ(0O) 
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Cent, 


1 


' u }<  M‘ 4* y (E+p • cosa-B2HC-  tana-p • sina)  (-4C)  (B^~S~^)  S-na  j 
0 [L  21  J 


K - M- 


2-^C 

21  J L 21 


+ (B2“E-C)  • sina  • (C* tana-p* sina) 


. < ' cl 

r4  J 

21  J 


*[  axl  *“B<B->  ' 1>w 


+ U2(B)  ' (MUB)  + U2  (B+H)  * (A4UPH) 


+ u 


2 (XI) 


b2+h2 


2*C* (K-M*  4*  u) 
(xT  - Bo)2 


(M*4*y(cos0-B2)2(l-B2*cos0) (-4C)  sin0/R6  ] 

23 


+ K_  B?)2  |-2*C*  sinS  + (1-B2cos0) 4C* sinfl  I ^ 

R2  * . L R R [ 


F 


■««(B)  ' <A*B>  * ' (MUP»> 


• / \ 


an 

V»2  (X1'B2>2 


3. 2. A. 4 


Itfrf  n 


(- 

■ I „5(p)  +M*4*m (E+p*cosa-  B2^(C" 

Jo  L 

V. 

[r  • . (E+o'cosa-  B?)^\  r 2 • C 

K-M-4 -P— V H ’ [ rZ” 

21  J L 21 

+ (B2-E-C)  • sina  • (C-tana-p-sina) 


.o.  ......  (B^-E-C)  sin  a 

(E+p-cosa-  B2r(C*tana-p-sinc-)  (-4C) — * 

21 


4 • C n 

“FJ 

21  J 


+ 1 v1(p)  <1  lM'4-y2*C(E+p*cosa-B2)  * 


• [p*cosa-(E+2C-B2)  +4  • (C*  tana-p  • sina)2 


(E+CHBo) 


+ K-  • sina  + [p+(E-B2) cosa-C-  tana  -sina] 

L 21  1 

• [c* tana-p *sinaj  (4C)/R^j| dp  +v2(b-H)  * ^A4VMH^ 


+ V2(B)  '(A4VB)  + V2(B+H)  ‘ (A4VPH) 


_ (cosG-  Bo) 

V3(0)  1M‘4^2*C 

[ 23 


jjl-B2’COs0+  sin20) 


-4  • (l-B2•cos0)• 


2C  2*B2Sinz0 


- cos0  ' d0 


23  23 


+V4(B-H)  * (A4VBH)  +V2(B)  ’ (A4VB)  +V4(B-H)  ' (A4VPH) 
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j V5(p) 


C(E+p*coso-B2)  - — 

21 


[?  (E+C-Bo ) 1 
p'cosa-(E+2C-B2)  + 4 . (C-tana-p *sina)  j 

K • 2 -2C  • sina  +[p+(E-B2) cosa-C*tana*sinaJ 

/■»  i ^ 


(c*tana-p*sinaj  (4C) /R^  j j 


dp  - X * FA4(B2) 


3.2 
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2*  B • [k-M’4*Ucos: 


+ u • (A5UMHP  + A5UMHN)  + U . * 


+ u (A5UPHP  + A5UPHP) 

1(B) 


+ 2 • B • [K-M*4*ucos2a] 


(A5UBP  + A5UBN) 


(1  (r 

+ J “2ml  m'4'"4' 

r [>- 


(Xl-E-B-cosa) 

R2 

12 


B * sina 


+ M‘4‘11  (Xl)*-E-B*cosa_)_2  (4.  sina.  [(E_3a+C)  B]  ) . C‘tan^B>  S1-n-^ 

R12  JL  12 


4 • sina 


+ u 


r (Xl-E-Bcosa)^J  2 •B’Cosa  _ (C:tana-B-sxnal 

■ Ik  - M-4*u  r2 gZ  r1* 

L 12  JL  12  12 

. 4 • sina  (E-xT+C)B)J  > dXl 

f Ifl-E-cose-C- tana -sin9-B» cos  (Q+a)]  . r4  . M.4.y~| 

+ Jq  U3(0)  jL  ^ 1 

r f f-nsfl-E-B  • cosa ) B • sina  . (cosG-E-B-cosafr-B- Sin(eH-a)+B^C),siiia)~| 

• h s2  r13  j 

L 13 

T m*4*u  (cos0-E-B’COsa) 2]  r -2-B*gin(9_+a) 

+ Ik j?  J|_  iT 

L 13  u 

(B» (E+C) sina-B* sin (9+a)fl  dQ 

- (l-E’cosG-C'tana-  sin0-B’ cos  (0+a)j*  4*  g1? J 
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r 


i 


«,®  • » • - 


+ M-4-MXl  -E-B-cosa)_2  (4<  s.na.  J-(e.^+C)bJ  >]  j C’ tanfl^- si-^] 
12  12 

T ( xi-E-B • co Sn ) 2]  f2‘B,cosni  _ (C^taiyi-B'sing) 

+ K-M"  4 * y ^2  ^ * R2  R4 

L 12  1 J L 12  12 

• 4 • sina (E-Xl+C) B)1 


dXl 


B~H~2 -b"k-M- 4 • u-cos2sJ  . d 

(p-B)2  5(p) 


+u5(b_h)  * (A5UMHP-A5UMHN)  + U5(Q)  • 
( A5UPHP- A5UPHN ) 


(A5UBP-A5UBN) 


+ U 


5(B+H) 


f -2-B-  Fk-M- 4* \|COs2a~[  # 

+J  (P'B)2 

B+H 


5(P) 


rB-H 

M*  4 • u • B • sin^a  . 

Vl(p)'  (p-B)2  ^ 


0 
+ v 


l(B-H) 


+ V 


1 (B+H) 


(A5VMHP  + A5VMHN)  + v1(0) 
(A5VPHP  + A5VPHN) 


(A5VBP  + A5VBN) 
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•j 


dp 


I 


v.  , ..  M*4*uB*sin2a 
K»>  ■ (p-B) ^ 


B+H 


* /\g a 


12 


-2-B-cosa  (-C- tana+B- sina)  (Xl-E-B-cosaJ/R4* 

12 


- (XI  -E-B-cosa)  (-C-  tana+B-  sina)2  -4-B  (E-Xl+C)  • sir*j/BS 

12 


r-B- sina 

(Xl-E-B-COSa)  -2-B-  (E-Xl+C)  • sirWl  1 

L r2 

R4  f 

L 12 

12  J j 

* J 2"v 3(e) 


' sina(sin0-C-  tana+B*  sina) 


(1-E-  cos0-C*  tana-  sin0-B*  cos  ( 0ta)  /R4 

13 


+B-cosa-  (cos0-E-B-cosa)  • (l-E-cos0-C*tana-sin0-B’cos(0+o)  )/R** 

13 


+ (cos0-E-B-cos«) (sin0-C- tana+B- sina)B- sin (©+0)/R4 

13 

+ [(cos0-E-B-cosa) (sin©-C- tana+B- sina) (l-E-cos0-C-tanasin0-B-cos(0+a) ) 


•4  • (-B • sin (&ta)+B(E+C) -sina)] /R6  j 

-2-K- |5_22£|®ia)  + (e- sin0-C-tana’cos0+B- sin(e+a) ) 

13  { 

-2- (-B-sin(€H-a)+B- (E+C) • sina)/R4  J> 


13 


d0 


3. 2. A. 5 Cont. 
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..  .1  « tnaA’iX  ™ ^ 


(-C*  tana+B-  sina)  2 


-2*B*cosa(-c*tana+B*sina)  (Xl-E-B-cosa)  /R*4 

12 

- (Xl-E-B’cosa)  (-Otana+B*  sina) 2*  4-B  (E-Xl+C)  • sina/R6 

12j 


+ 2-K 


-B*  si 

I i*2 


sina  (Xl-E-B-cosa) •2-B- (E-Xl+C) -sina 


+ V, 


(-M«4-pB‘ sin2a) 


+ V5(B_H)  * (A5VMHP  - A5VMHN)  + vg  • (A5VBP  - A5VBN) 


+ V5(B+H)  ' (A5VPHP  - A5VPHN) 


(-M-4-vi-B-sin2a)  , 1 

* v5(pl ap  I ‘ FA5(B) 
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( B-H 


u 


lip)  P-B 


<=?>  ‘<>  * “KB-H,  • (B1UMH)  + “1(B)  • (E1UB> 


+ U v * (B1UPH) 

U1(B+H) 


/'C 


,+K  , , + u — ~ Jm-4'u^-E-B-coso] 

ui(P)  ‘75’  p J 21x11  I 

B+H  F 

[B-sina-C-  tangl  2_  R . [xL~E~^‘  c9..?.°lI  j dxl 

R4  19-1 

12  12 

+ f u3{0)  |-M-4*y[cos0-E-B.cosa].[sine-C-tana+B'sina] 

0 '■ 

• [1-E-  COS0-C-  tana-  sin0-B-  cos  (O+a)]/^  ^ 


-K-  [E-sinQ-C 


- tana- cos0+B»  sin (0+a>]  /R*  o jd0 


. f1  u — i -M*4"  p fxl-E-B'COsaJ'jB'  sina-C- tsinot]  2/r' 
4 (Xl)  L 


12 


+ K • [xl-E-B-cosaJ/R^|  dxl 


B-H 


1 U5(p)  ‘F5>  d>  -“mb-h,  • (B10“H  - "mb,  • 


K 


(BlUB) 


"u5 (B+H)  • (B1UPH>  + J U5(P>  'P“B 

B+H 


dP 
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126 


f 


) 

l(p)  1 


■M*  4"p  [EMp  •cosa-B2]-[C'tana-p  ' sina] 


• f (E+p  *coso-B, ) (-sina ) + (C- tana -p  "sina)  (-cosa)J  /Bf* 
u 21 

-K  • |e+p • cosa-B2 ) (cosa)+  (C* tana -p  • sina)  (-sina)]/R2  | dp 

rB2-«2 

J u2  (xi; 


-K 

:(xi)  vxi-B2 


{^7  > dXl  + u2  (B-H)  * (B2UMH)  + U2  (B)  * (B2UB) 


+ U2 (B+H,  • (B2°PH)  + 


j u2(xi: 


{ 


(XI)  1 X1-B2 


} dXl 


B2+H2 


f2”  I 

I -M-4'u  (COS0-B,)  (sine)  (I-B-cosqJ/R1* 

J 3(0)  | H 2 23 


-K[B2-sin0]  /R2^|d0 


B2-«2 


+ | U4(W)  {lfc 


B2}  dXl  “ U4(B-H,  * (B2UMH)  - U4(B)  ' (B2°B) 


' U4 (B+H)  ' (B2UPH)  + 


r • 


4(3  axl 


V»2 


* (\ 


5(p) 


-M’ 4*  p [E+p • cosa-B2]-[C’ tana-p  * sina] 


+ [(E+p-cosa-B2) (+sina)  + (C-tana-p-sina) (+cosa)]  /R^ 

1 


+ K [(E+p*cosa_B2)cosa  +(C* tana-p’ sina) (-sinai^R2^  | dp 


3.2.B.2 


i 

I 


I 
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IM*-  + |o  |>-|>- 


+ v. 


(C- tana-p • sina) 2 


(jE+p* cosa-B2) (-sina)  + (C* tana-p* sina) (-cosa)l/R2  Idp 

21 


1 

— • v + 

X_  2(B2) 


f 2lT  „ (r.  M*4*p  (sin26)l  f l-B?cosQ  1 1 

j 3,e)  J'L  ^ 'J j 


4<B2) 


v5(P)  r mm‘11 


(C* tana-p* sina)  21 


>*cosa-B2)sina  + (c* tana-p* sina )cosa]/R2  idp 

21 
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• [k*sina+c*sina-sin  (a+0p)J/R^|  dp 


I*  FB3(0O) 
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/ V.,  ![■ 


• 4‘  ]i2 • C (E+p  • cosa-B2) 


sina 


21 


• [p*cosa-(E+2*C-B2)  + 4*  (C’tana-p*  sina)  • 1 

21  i 

-K  *sina+  [p+(E-B2)  •cosa-C’tana*  sinaj 

L 21 

• [c* tana-p • sina]  (4C)/R4  J dp 

+ U2  (B-H)  * (M'JMB’  + “2(B)-  ,BWB)  * U2  (B+H)  * <MUP“> 
( 2lT  / 

hj  U3(Q)  |m*4*U*  (2C)  • j\-B2cos0+sin20) 


sin20  1 „ 2-C  T 2*B9*sin20  „"|  i 

s'-  J ’ L ~ P COS0J  [ 39 

23  23  23  J 


,,,  sin 
-4*  (I-B2  * cosO) — -2 
R 


+ U4(B-H)*  (B4UMH)  + U4(B)’  (B4UB)  + U4(B+H)’  (B4UPH) 


5 (p) 


M*4*p2*C(E+p*cosa-B2) 


21 


2 | B 0 ) "1 1 

• [p*cosa- (E+2«C-B2)+4* (C* tana-p* sina)  • R2  7 j 

21  J 

-K  *sina  + [p+(E-B2) *cosa-C*tana-sina3 


[c* tana-p* sina]  (4*C)/R4 
u J 21 


dp 


3.2.B.4 
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.. 


j-  : -v-.a,  . .. 


a 

i 


/ vi to)  ![&•«•»• 


(B2-E-C)'sina*  (C*tana-p- sinct)  (4c)J  • ^ 


ft  _ (C*  tana-p • sina) 2 
L R2" 

21 


+ tana-p » sina) 2 -| 

21 

’ [p+  (B2-E-C)-sino-  (C-tana-p.sinaJ-^-JjJ 


/'  b2~h2 

+ / v — • -(+K.2.C)  — 

J 2 (XI)  <Xl-B2)2  dX1  + V2  (B-H) * (B4VMH)  + v2 (B) * (B4VB) 


+V2 (B+H)  -<b4VPH)  + 


v — . 1+K‘2'C) 

2 (XI)  (Xi-B2)2  dX1 


b2+h2 


2, 

+ V„ 


, 2 • C* sin0 


M-4-p.l^in0  . 

R4 

23 


T,  sin20  1 

I ‘ 

1-B2*cos0| 

23  J 

L J 

23 

B2_H2 


23  23  J 


+ V — 

y 4(xi: 


(~K*  2*C)  , — 


(XI)  (Xl-B2)' 


dXl  - 


^4 (B-H) * (B4VMH> 


~ v4(b)  • (B4VB)  - v4(b+h).  (B4VPH) 
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B-H 


u,  M* sin2a 

l(P)  (-D-B) 2 dP 


+ U1  (B_H^ • (B5UMHP  + B5UMHN)  + • (B5UBP  + B5UBN) 


+ U1(B+H)  ‘ (B5UPHP  + B5UPHN) 


'/  *> 

B+H 

f 

I „ ..  - T_  (-C*  tana+B*  sina) 2 

■ U2  (XI)  l 2*M,4,y  |-B*sxna ^ 

J l 12 


-2*B,cosa(-C*tana+B*sina) (Xl-E-B*cosa)  /R4 

12 

- (Xl-E-B*  cosa) (-C* tana+B*  sina) 24 • B (E-Xl+C) • sina/R6 

12  -* 

_ ..  r~B*  sina  (xT-E-B’cosa)  2*  B (E-Xl+C)  sinal  (.  — 

~2*K  L ? j dxl 

12  12  -'J 

+ [ 27F  f r 

J u3{0)  j 2•M•4•^l•  |b* sina (sinQ-C* tana+B* sina) 


• (1-E‘cosO-C*  tana*  sinG-B*cos  (0+a)/R4 

13 

+ B*cosa (cosO-E-B'cosa)  (l-E'cosQ-C'tana^sinG-B’cosfG+a)  )/R4 

13 

+ (cosO-E-B*cosa) (sinG-C* tana+B*  sina) B* sin (0+a)/R4 

13 

+ [(cosG-E-B*  cosa) (sinQ-C*  tana+B*  sina) (1-E*  cosO-C*  tana*  sinQ-B*  cos (0+a) ) 


3.2.B.5 


* 4(-B*sin(0+g)+B(E+C)sing) 


+2*k  ~~  ^~~a^ ■ + (E* sin0~C,tana*cos0_B* sin  (0+ot) ) 

R 


*2* (-B*sin(0+g)  + B* (E+C) sing)/R4  id© 


■ J U4(3a)  !+2*M‘4*»1[-B*. 


(~C»tang+B«sing)  2 


-2*B*cosg(-C*tang+B*sing)  (Xl-E-B*cosg)/R4 

12 

- (Xl-E-B*cosa) (-C*tang+B*sing) 2*  4*B (E-Xl+C) •sina/R6 J 

„ „ f-B*sina  (Xl-E-B*cosg)  • 2 *B  (E-Xl+C) sing! 

-2-K  1—^2 & -I  dXl 

12  12 


(-M»4,u*B*sin2g) 


- dp  + u_.  . • (B5UMHP  - B5UMHN) 

5 (p)  (p-B) * K 5(B-H) 


+ , . • (B5UBP  - B5UBN)  + Uc . • (B5UPHP  - B5UPHN) 

5(B)  5 (B+H) 


(-M»4«u«B»sin2g) 


r°-"  ] ] 

+ J Vl(p)  j 2 * B*  4 * P • sinzgj  ~J^7  j 


+ V,  ...  • (B5VMHP  + B5VMHN)  + V,  • (B5VBP  + B5VBN) 
l(B-H)  1(B) 


+ v,  . • (B5VPHP  + B5VPHN) 

1 (B+H) 
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J Vl(p)  i2*B*[K-M*4*p*sin2a]  -fiZSp  <iP 
B+H  *■  J 

J V2(xT)  2*M*4*y*J(B*sina-C*tana)  *2*B*c 
F l 


+ (B*  sina-C* tana) 2 * 2:Bvsim  (^C-Xl)j J"  C-_t_ana-B* sina 

12  J L 12 
(B* sina-C* tana) 2 


2*B*cosa  . 2*  (B* sina-C* tana) *2* B* sina (E+C-Xl) , , 

r? + rn f dxl 


+ v. 


I r 

|2*M*4*p  |^( sin0-C* tana+B* sina)  *2*B*c 


. r . ~ _ . i?  , . [(E+C)  sina-sin(a+0)H 

+ |_sin0-C*  tana+B*  sxnaj  z*  2*  B*  v ^7-  -- 1 

13  J 

• 1-E* cos0-C* tana* sin0-B* cos (0+a)  /R4 
L J 13 

. T.  ..  „ (sin0-C* tana+B* sina) 2 1 

+ K - M*4*y  £2 

L 1 


2— ^ siri(Q+a) — 2*  (1-E* cos0-C* tai 
13 

•2*B*  ( (E+C)  * sina-sin(a+0)  )/Rl*  ' 

13  J 

^1  ( J 


2* (1-E* cos0-C* tana*  sin0-B*cos(0+a) ) 


V4  (XI)  i+2*M*4*y*  I (B*  sina-C*  tana)  • 2*  B*  cosa 


+ (B* sina-C* tana) 


2 *2*  B*  sina (E+C-Xl)] | C*  tana-B 


1 C* tana-B* sina  1 

f J 
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f 


(B*  sinct-Otanct) 2 


[2*B*cosa  2*  {B* sina-C*tano)  • 2* B* sina (E+C-Xl)l 

i2 + ^ j 

12  12  J 

+ J v5(p)  |-2.B[K-M-4-ysin2a]*-^Ii)2|  dp 


+ v_,„  • (B5VMHP  - B5VMHN)  + v_/fa  • (B5VBP  - B5VBN) 


+ v5(B+H)*  (B5VPHP  - B5VPHN) 


-2*B*  [k-M*4*n*sin2o]*  1 dp 


I * FB5(B) 


3.2.B.5  Cont. 
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Appendix  B.2 

Algebraic  Insert  Formula 


B.  2 Insert  Formulas  for  Eq.  set  (3.2) 
Insert  Formula  for  Linear  Displacement  Assumption 

(u  = p + c2) 


AlVMH  5 -K 
AlVB  E 0 

AlVPH  E +K 
Linear  Insert  Eq.  A2 

A2VMH  E 

-K 

A2VB  2 

0 

A2VPH  E 

K 

Linear  Insert  Eq.  A4 

A4UMH  2 

0 

A4UB  5 

-4  *C  • (K 

A4UPH  S 

0 

A4VMH  E 

it  • C • (M 

A4VB  E 

0 

A4VPH  E 

i 

=i 

• 

o 

• 

Linear  Insert  Eq.  AS 

A5UMHP  E 

M*  2 • y •: 

A5UBN  E 

16 ‘M*  n 

A5UPHP  E 

-M  • 2 • y 

A5UMHN  E 

A5UBP  * 

A5VMHP  2 

(M»2 

A5VBN  2 

-M»  8 • n 

B 

'h' 


B 

’h 


TT_ 

’h 


Bit 

► 

H 


B 

> — • 
H 
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.***.*1  ■ iaMr^tiuVr-r  


r-~  -i  » 


A5VPHP 


A5VMHN 


- (M*  2 »)i  • cos  2a  - K)  ♦- 
A5VBP  = A5VPHN  = 0 


Linear  Insert  Eq. 


BlUMH 


B1UB 


BlUPH 


Linear  Insert  Ee 
B2UMH  = 


B2UB 


B2uPH 


Linear  Insert  Eq. 


B4UMH 


TT*C*  (M*2  + K)/H, 


B4UB 


B4UPH 


-n*C* (M*2u  + K)/H, 


B4VMH 


B4VB 


B4VPH 


Linear  Insert  Eq. 


-4  *C*K/H2 


B5UMHP 


(M*  2 • v 'cos  2a  + K) 


B5UBN 


B5UPHP 


-M»  8 •==■•  sin  2a 

-(M*2*ycos  2a  +K)  • 


B5UMHN 


B5UBP  = B5UPHN  = 0 


B5VMHP 


B5VBN 


B5VPHP 


-B*M*2*y*g-  *sin  2a 
4*  g-*  (M*  4 • y • sin^a  -K) 
B'M^^y*  j7*sin  2a 


i 


B5VMHN 


B5VBP  ~ B5VPHN  = 0 


Insert  Formula  for  Square  Root  Displacement  Assumption 

(u  = c.  1/2  - ' 


+ c2) 


Define : 


n/b  + If  - 2\fB  +\IW=TT 
(JB+H  + \fB ) • (\/B-H  - n/B) 


Square  Root  Insert  Eq.  A1 


A1VMH  E k\2\[B  - 2n/5=R  - . 5^  *Ln  (AG)  ] D^ 

AlVB  E K [2'v/B-H  - 2\[B+H  + n/B -Ln  (AG)]  ~ 

A1VPH  E K |-2\/B  + 2\lB+H  - . 5^»Ln(AG)|  1 

D1 

Square  Root  Insert  Eq.  B1 


AlVB 


A1VPH 


BlUMH 

B1UB 


BlUPH 


- A1VMH 


-AlVB 


- A1VPH 


Square  Root  Insert  Eq.  A5 


A5UMHN 


A5UBN 


A5UPHN 


A5UMHP 


A5UBP 


A5UPHP 


A5VMHN 


A5VBN 


(M*4  *vi  *cos2a-K)  Ln(AG)  + —) 

(M*4  *y  «cos2a-K)  (-  *jB  Ln(AG)) 

A5UMHN 

-M*  2y  ^l_*sin  2a  (H  + 2B  - 2'jB'M+R) 

n*  U* 

M A "TBIT  

-M«4*yH<Di>sin  2ct  (-H  + «s/B\/ I+fi  - s/Bn/iFR) 


-M*  2 • 


-sin  2a  (H-2B  + 2*jB\IW=il) 


H R Wr 

-M*  4 •sin12a  (g-  + — ^■•Ln(AG)) 

Jo 

M*4«p*sin  2a  »Ln(AG)) 


I 

! 

i 

R 

! 


A5VPHN 

- 

A5VMHN 

A5VMHP 

— 

IT 

(K  - M*2  *y  *cos  2a)  (2  •B/B+H/B:-2BV/B+H 

H *D^ 

A5VBP 

TT 

(K  - M*2*ycos  2a)  (2*  B/B+H 

H*D, 

1 

- 2/b»H  - 2 *B  /B-H) 

A5VPHP 

= 

TT 

(K  - M*2«ycos  2a)  (2  •b/b-H 

H #D, 

1 

-2*Bf§+  H/B) 

Square  Root 

Insert  Eq. 

B5 

B5UMHN 

= 

A5UMHN 

B5UBN 

A5UBN 

B5UPHN 

= 

A5UPHN 

B5UMHP 

~~ 

-TT 

(K  4 M*2*ycos  2a) 

H*D1 

(2*B/B  - 2-BjB+H  + hJb) 

B5UBP 

-TT 

(K  + M*2*ycos  2a) 

H *D, 

1 

(2*b/b+h  - 2/B*H  - 2 •b/"b-h) 

B5UPHP 

- 

- IT 

H*D, 

(K  + M»2*y»cos  2a) 

1 

(2  •B/B-H  - 2*b/b  + H/B) 

B5VMHN 

(M*  4 • 

y • sin2a-K) (/b  *Ln(AG)  + 2B) 

20^^  H 

B5VBN 

(M*  4 « 

y*sin^a-K)  (-/§*Ln(AG) 

B5VPHN 

D1 

B5VMHN 

B5VMHP 

M*  2y  • 

2a  (H  + 2*B-  2 /b/b+H) 

B5VBP 

5 

M*  4y  • 

——•sin  2a(-H  + JbJb+H  = /E/b-H) 
H-D1 

B5VPHP 

K*  2 • y 

•JB'V  sin  2a  (H  - 2*B+  2JbJb=H) 

I 


Insert  Formula  for  General  Non-Linear  Displacement 

Assumption 

(U  = M I-  C^} 

Solution  of  the  following  transcendental  equation 

sin  y Ctt  + a] ) +y  sin  [tt  + a]  =0 
Define: 

d2  [C  - B + H]y-  [C  - B - H]y 

ill  -1) <C  ~B) Y~2 

2 *D- 


D21  ~Y(C  -B)  y~1  — B • Y (Y  -1)  (C-B)  Y~2 


— 1 

Y-  Insert 

Eq. 

A1 

D2 

A1VMH 

= 

”D11*K*  4 *B*H  “ D21*K*2*H 

A1VB 

= 

0 

A1VPH 

= 

-A1VMH 

Y - Insert 

Eq. 

B1 

B1VMH 

= 

-AlVMH 

B1VB 

= 

0 

B1VPH 

= 

-AlVPH 

y - Insert 

Eq. 

A5 

A5UMHP 

= 

2 

”D11  ( B *M*  ^ * V [4  ’cos  a,H-2*B*,rr,sin2a] 

-K*B*4*H) 

A5UBP 

B 2 B 

-M* 4 •ir,cos  a+  K’jtM 

n n 
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A5UPHP 

= 

-A5UMHP 

A5UMHN 

= 

A5UBN  = A5UPHN  = 0 

A5VMHP 

- 

D21*  (M*4  *B,7T  *cos  2ct"2*K*7I*B) 
+D11*(M,8*M*B(H*sin  2a+  b*tt*cos  2ct) 

-2*tt*b2*K) 

A5VBP 

= 

M*  8 * U *77  sin  2a 

n 

A5VPHP 

= 

-A5VMHP 

A5VMHN 

= 

A5VBN  = A5VPHN  = 0 

Insert 

Eq . 

B5 

B5UMHP 

= 

D2i*  {M*4*li*B*Tt*cos  2a+  2*K*n*B) 

+ D^^*(M*8*M*B(H*sin  2a+  B*n*cos  2a ) 

+ 2 * ^ *B2  * K ) 

B5UBP 

= 

M*8*y|*sin  2a 

B5UPHP 

= 

-B5UMHP 

B5UMHN 

B5UBN  = B5UPHN  = 0 

B5VMHP 

= 

-D11*(B*M,4,U(4,sin2a*H  + 2*B*Ti,sin2a 

-K*B*4  *H  ) 

-D21*M'4  ’TT'S*  sin  2a 

B5VBP 

= 

B 2 K 

-M*  16*  V ’{t*  sin  a+  4 *B*- 

H H 

B5VPHP 

= 

-B5VMHP 

B5VMHN 

B5VBN  = B5VPHN  = 0 

Define : 

°3  tB2  + H2  _F]Y  ‘ tB2  " H2  "F]Y 

C11  ■Y(Y_1) {B2  ‘ F)Y~2 
2D^ 

C21  _Y(B2  “ F)Y"1  + b2*y*(  -1MB2-F)y_2 
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r 


y-Insert  Eg.  A2 


A2BMH  = 

C11‘K*4*B2  V C21*K-2*H2 

A2VB  = 

0 

A2VPH  = 

-A2VMH 

Y-Insert  Eq.  B2 

B2UMH  = 

-A2VMH 

B2UB  = 

0 

B2UPH  = 

-B2UMH 

y-Insert  Eq.  A4 

A4UMH  = 

+Cxl • 4 *C (K  - 4*M*y)H2 

A4UB  = 

-4 *C (K  - M*4*y)/H2 

A4UPH  = 

-A4UMH 

A4VMH  = 

(Ci;l*2.B2  + C21)2*tt*C(K  - 2*M»y) 

A4VB  = 

0 

A4VPH  = 

-A4VMH 

Y-Insert  Eq.  B4 

B4UMH  = 

-(C11*2*B2  + C21)2*tt-C(K  + 2*yM) 

B4UB  = 

0 

B4UPH  = 

-B4UHM 

B4VMH  = 

C11*K*4*C*H2 

B4VB  = 

-K*4*C/H2 

B4VPH  = 

-B4VMH 

E 


I 

t 
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■>.  f-  S-  ■ i.M.  ■■  -■■„  . J. 


f 


( C 


B-H 


(txl(p)  _txl  (B)  * 2 •Ln<P"B)2dP  + fcxi  (b)‘  2 


Ln (p-B) 2dp 


+ Scl  (B)  [2-H‘LnH-2.H]  + I Ln(p-B)2 


xl (B)  2 


B+H 


+ J ^xKp)  -M-cosZta)  - tvl^rM.sina*cosa]  dp 


yKp) 


♦/  (tx2(xT)'  2 ' tn(RJ2)  + * 


x2(Xl) 


M. (Xl-E-B.coso) 2/R2 
12 


+ t — • M- (Xl-E-B-cosa)  (B»sina-C«tanct)/R2  } dXl 
y z ( X J. ) 1 2 

2* 

+ | {t  . • ^ • Ln*R2  + t . • M* (cose-E-B«coso) 2/r2 
I xJ  2.  23  x3 10;  2 


+ t *M* (cos0-E-B*cosa) (sin0-C*tana+B*sina)/R2  } d0 

13 


+ {tx4(Mj  ‘ 2 ’ ■ «■  W-E-B-cw,)*/* 


+ t — .M*(Xl-E-B.cosa) (B.sina-C.tana)/R2  } dXl 
i*»i ) 22 

, C , B-H 

+J  (tx5(p)"tx5(B))’ 2 ,Ln(p_B)2dp  + ScSCB)'!  J Ln  (p-B) 2 dp 


+ fcx5(B)  i2-*-™-2-*! + 1 


x5(B)  2 


Ln(p-B)2  dp 


B+H 


FA1 (B) 


+ tx5(p)  • M*  cos2a  - ty5(p)  .M  .sina*  cosaj  dp 
J 0 
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- .x,  n-a.li  --  x.  _ .Uj*.  - 


-4-^- 


I (txl(p)'  2 • R2,  * ‘xl  (p)  * "•  (Wp.co.rBjlVn^ 

0 

+ tyl(p)*  M • (E+p.cosa-B2) (C.tana-p.sina)/R2  } dp 


+ J { (tx2  6a)  -fcx2  (B2) 5 • I * (X1-B2>  2 + Sc2  («L)  • M > dX1 


+ fcx2(B2)  J J - Ln(Xl-B2)2  dXl  + tx2(B  2 (H.Ln  H-H) 
F 2 


+ t 1 — 2 — 

x2(B2) J |-Ln(Xl-B2)  dXl 


. 2ir 

+ J {tx3(0) ’ f Ln  R23  + fcx3 (q) * M.(COS0-B2)2/R23 


+ t • M.  (coso-Bp)  (sino)/R2  } do 

yj  <•  o q 


* (<tx4<5T)  -W)'-  |-i««i-b2)2  ♦ tx4(-r  H,  axi 


1 2 

+ — • Ln  (Xl-Bo ) dXl  + t _,D  2.  (H.  LnH-H) 

x4(B2)  I 2 2 x4(B2) 

; F 


+ tx4(B2)f  i-Ln(X1-B2>2  M1 


+ I {tc  • ^ *Ln  R2  + t *M  (E+p.COSa-B2)2/R2 


x5  2 21  x5 


21  FA2(B2) 


+ t _•  M •(E+p*cosa-B2)  (C*tana-p»sina)/R2  } dp 

y=  2i 
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f 


tScKp)  * 2 *ln  R^'+  txl(p)*M-(E+P*cosa-co8e0)2/R2i 


+ M* (E+p*cosa-cos0o) (C.tema-p.sina-sin0o)/R2  } dp 


*r 


r “ R32  + Cx2 (xT, * 


t _ • (Xl-cos0o) (-sin0o)/R2  } dXl 
y2  u 32 


f2* 

+ j (tx3(0)‘  "•Si”2  * 'y!(ei  • I • si”(e*B0>) 


2tt 


+ I < x3<9)-2  x3<.9o}  . Ln(2-2.cos(0-eo)) 

0 

/-  0o_Ae 


d0  + 


^(Onf  2 


Ln(2-2*cos(0-0o)d0  + (q  ^ • I Ln(2-2.cos(0-0o))d0 


■»■/ 


0O+A0 


+ t 


A0 

x3(0n)  * 2 


. [4.Ln(M))-2.Ln.4-4  - -^-2  ] 


+ f {^2,  /~x  * T,Ln  rZ  + * M>  <X1,  o2SG'n) 


x4 (XI)  2 32  x4(Xl) 


32 


+ t . — • M • (Xl-cos0o) (-sin0o)/R2  } dXl 

y4\Al;  v 32 


+ f {tx5(p)‘i,Ln  R31  + ^(p)*  M*  <E+P-cosa-cos0o)/R23i 
' 0 


+ t • M*(E+P’cosa-cos0o) (C*tana-p*sina-sin0o)/R2  ) dp 
y5(p)  w 31 

FA3 
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t 

I 


t . . . (-2  C) (C*tana-p*sino)/R2 

x±ip;  21 


+txl(p)*  M ‘ <E+P  *cosa-B2)  2 (+4C)  (Otana-p^sinoO/R1* 

M.^(E+p.cosa-B2)]^|^-  + ,(CI.tana-p|>sina)2(+4C)j 


+ t 


yi(p) 


dp 


( B2-H 

r -M*2-C  ] 

J y2  (XI) 

L xi-b2J 

F 

ftx2(B2>  &'C 

(M-l)  • TrJ 

( 2"( 

I.  r- 

■2*C*sin0  1 

J | X3 (0)  [_ 

"TP J 

23 

B,+H 


-B2 


dXl 


23 


r -2’C 

sin20  (4C)  1 

R* 

23 

23 

d0 


x4  (B2) 


(j2 * C (M-l)*ir] 


,b2-h 

J 4-  ■ i - 

r-M-2-c~| 

dXl  + 

fl  t - 

-M*2*cl 

J y4(xi) 

L xi-b2  J 

/ y4 (XI) 

521 -b2J 

dXl 


B2+H 


x5  (p) 


-2«C«(C,tana-p,sina) 


21 


FA4(B2) 


LJ 


151 


' fcxl(B)  ,2‘B  * ff+  fcxl (p) * M‘sin2a  dP 


B-H,C 


t • M » sin2a»- 

xl(p)  p-B 


- • M • 2 • B »ir»  sin2a 


0,  B+H 


+ { t • M • cos2a  dp 

yi(p)  w 


B-H,C 


t • M • eos2a,-^-fB-  dp 

yl(p)  p-B  M 


0 jB+H 


“ tyl (B)  *M  • 2 • B • it*  sin2a 


+ t — J»2-B*sina  (C+E^X1)  + M-m-ErB^cosa)  2(X1-E.-C).4  • B-  sine 

I X*  v XX ) I R 

JF  l 12  12 


-4*M-(+B*sina) (Xl-E-B-cosa)/R2  IdXl 

12 

. ^ (B*  sina-C* tana) . ~ w(Xl-E-B*cosoO  , _ 

I fcy2 (xl)  i“M*2'B,sln“  £2 h2*M 52 (-B-cosa) 

■'f  { 12  12 

+M  • (Xl-E-B*  cosa)  (B • sina-C* tana)  4— B Sin°4  ,E  1 dXl 


( 2ir  ( 

+ t < 2 - D ^ fE+cl  sina-sin (a+0) ) 

x3 (0)  I B r2 

J0  y 13 


-4»b*M»  sina (cos0-E-B*cosa)/R2 

13 

+ M*  (cos0-E-B* cosa) 2»  4 • B*  {sin (a+0) - (E+C) sina) /R1*  l d0 

13. 


+ I \ -2»M»b* sina*  {sin0-C*tana+B»  sina)/R2 

J y3'0'  13 

0 { 


FA5(B) 


■ ,4* 

I 1 
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2»M*  (cos0-E-B* cosa)* B* cosa  /R2 


13 


+4*B*M* (cos0-E-B* cosa) (sin0-C* tana+B* sina) 
1 


. (sin(a+Q)~ (E+C)  sina) 


inal  1 


13 


f 


1 


x4(Xl) 


«.B-.ln.(C+E-«l)  . (4.B.sim)  ma^S2^l 

R , * 

12  12 

] 


+ M (Xl-E-B'cosa) 2*  4*B*  sina (Xl-E-C)/R^ | dXl 


+ 1 t . I 2«M* 
J y4  (XI)  | 


. (B*  sina-C*tana) 
(“B*  sina) £2 


12 


+ 2*M* (Xl-E-B*cosa) (-B*cosa)/R2 

12  — 

(Xl-E-C) 

+ M (Xl-E-B* cosa) (B* sina-C* tana) (4 "B* sina)  ^ 

12 

,B-H,C  ’ 


axi 


+t  c *2  • B * ir*  + 
x5(B) 


JQ , B+H 


t . *M(-1)* sin2a  (p+p}  dp 
x5(p)  p-B 


t e/  • M* sin2a  dp  + t •**  • 2-B*7r-sin2a 

x5(p)  xd^b; 


f B-H.C  j 

t , „ 

-M*cos2a  i 

J y5(p) 

p-B 

0 ,B+H 

* J 

dp  + Vcb)  EM*2,itB,sin2aJ 


♦ f"  t 


y5(p) 


M (cos2a  ) f dp 


FA5 (B)  Cont 


A 


B-H,C 


V(P)'  2 -I,n(P"B,!  dp  + tyl(B)  *2*  (H’LnH-H) 


0,B+H 


+ Vtp)  *M‘sin2(a) 
} 0 


dp  + 


ScHp)  -M-(-isin2a)  dp 


j |ty2(Xl)*  2 


2 + t _ . *M*  (B*  sina-Otana)  2/R2 

12  y2(Xl)  12 


+ t , *M  • (Xl-E-B* cosa) (B*  sina-C*tana)/R2  > dXl 
(XI)  12 


+ f ft  R2  + t , . *M*  (sin0-C*tana+B*  sina)  2/'R2 

I y3  (0)  2 i3  y3(0)  13 

J0  { 

+ t _ .*M • (cosO-E-B’cosa) (sin0-C'tana+B* sina)/R2  1 d0 

x3(Q)  13 


♦ (l!t  - -A. 

j ] y4(Xl)  2 


— >Ln  R2  + t . . ’M*  (B*  sina-C'tana)  2/R2 
2 12  y4(Xl)  12 


+ t . ; *M» (Xl-E-B* cosa) (B* sina-C*tana)/R2  l dXl 
X4 (XI)  12  f 


, B-H,C 

+ t • Ln (p-B) 2 dp  + t •2»(H»LnH-H) 

I y5(p)  2 y5(B) 

^ 0 , B+H 


. C C 

+ t _ . .♦M*sin2a  dp+  | t c.  •M*(-  sin2a)  dp 
1 y5(p)  x5(p)  2 

J0  ' 0 


FBI (B) 
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I ~k  .1 

L i yl((”  2 


•Ln  R2  + t „ , *M* (C*tana-p* sina) 2/R2 

21  yKp)  21 


+ t . *M  • (E+p • cosa-Bo ) (C*  tana-p • sina) /R2  > dp 
xl(p)  z 21 


B-H,l 


',!|i  •yLn(X1-B2)2  <3X1  + tyJ(B2).2.(H.LnH-H) 


F,B+H 


,2ir 


+(  I fcy3  (0)  *2 


-•Ln  + t 


23  y3(0) 


•M  • (sin20)/R2 


23 


+ t *M* (cos0-B?) (sin0)/R2  > d0 

x3  (0)  z 23 


- B-H, 1 


Vmi  T“«i-B2)2  m + V(B2)'2'<H’LnH'H) 


F,  B+H 


+ It  R2  + t ...  , *M  *(C* tana-p*  sina) 2/R2 

y5(p)  2 21  yS(P>  21 

r0 


+ t _.  . *M*  (E+p* cosa-B2) (C* tana-p* sina )/Rz  l dp 

x5(p)  z 21  I 


FB2 (B  2 ) 
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0 


jvtp)*  2*Ln  R3X+  tn(p)‘M‘(C'tanC,"P‘SinCl"Sin0o)2/R31 


+ fcX2 (p) (E+p-cosa_cos0o)  (C-tanct-p-sinct-sin0o)/R2  ’dP 


' L / 
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figure  8.  Reference  Point  Locations  for  Bent  Crack  Geometry. 


Figure  9.  Pressure  Loaded  Crack. 
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Figure  I Ob. 


Figure  10c.  Traction  Loading  on  Bent  Crack. 
Figure  10.  Principle  of  Linear  Superposition. 


